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1 Introduction 


The evaluation of multi-loop Feynman diagrams has a long history [|^ and, roughly speaking, 
we can say that there are theories more simple than others, noticeably QED (see for instance 
ref. [0) and also to some extent QCD (see for instance ref. |Q), where we have few masses and 
the analytical approach can be pushed very far. Conversely the full electroweak Lagrangian shows 
several masses, ranging over a wide interval of values, therefore making the analytical evaluation of 
Feynman diagrams a complicated task. Equivalently complicated is the situation in any extension 
of the minimal standard model. 

Soon or later the analytical approach will collapse and one can easily foresee that this failure 
will show up at the level of a complete two-loop calculation in the standard model or beyond. This 
calculation is, for instance, required to produce quantities as sin^ [|| with a theoretical precision 
of 1 X 10“®. For this reason one is lead to consider an alternative approach to the whole problem, 
namely to abandon the analytical way in favor of a fully automatic, numerical evaluation of 
multi-loop diagrams. The Holy Graal (General Recursive Applicative and Algorithmic Language) 
requires fast and accurate procedures to deal with the singularities of an arbitrary diagram. 

For a fast and accurate numerical evaluation of arbitrary multi-loop Feynman diagrams there 
is some interesting proposal in the literature 0 that has been recently applied to the two-loop 
two-point sunset topology 0. Older approaches can be found in ref. p. 

The Bernstein-Tkachov theorem (hereafter BT) tells us that for any finite set of polynomials 
Vi(x), where x = (xi,... ,Xjv) is a vector of Feynman parameters, there exists an identity of the 
following form: 

V (x,d) n vt‘*\x )=B n vnx). (1) 

i i 


where P is a polynomial of x and di = djdxi] B and all coefficients of V are polynomials of /ij and 
of the coefficients of Vi{x). Any multi-loop Feynman diagram can be cast into the form of Eq.(|l|). 
Iterative applications of the BT functional relations, followed by integration by parts, allows us 
to express the integrand in parametric space as a combination of (polynomials)^ x logarithms of 
the same polynomials, therefore achieving a result that is well suited for numerical integration. 
The B coefficients of Eq. (|I|) will contain all Landau singularities 0 of the corresponding diagram. 

For general one-loop diagrams we have an explicit solution for the polynomial V which is due to 
F. V. Tkachov (see also ref. and ref. |^). Any one-loop Feynman diagram G, irrespectively 
from the number of external legs, is expressible as 


G = [ dxV-^{x), 

Js 

where the integration region is Xi > 0, < 1 and where V(x) is a quadratic form of x, 

V(x) = X* H X + 2 K'" X + L. 


( 2 ) 

( 3 ) 


The solution to the problem of determining the polynomial V is as follows: 

V = 1- X = K^H-\ B = L- K. 

2 (/i+1) ’ 


( 4 ) 
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For iV-loop diagrams with N > 2 a minimal BT approach has been recently proposed in ref. |3 
where we have adopted a different strategy aimed to deal with arbitrary diagrams. It represents a 
compromise based on the simple observation that we know how to apply the BT-iterative procedure 
for arbitrary one-loop diagrams. Therefore, given any two-loop diagram G we apply the BT 
functional relation to the one-loop sub-diagram of G which has the largest number of internal 
lines. In this way the integrand for G in x-space can be made smooth, a part from the factor B of 
Eq.(|ID which is now a polynomial in Xg, the Feynman parameters needed for the complementary 
sub-diagram of G with the smallest number of internal lines, Gg- The sub-diagram Gg, after 
integration over its momentum, becomes an additional - Xg-dependent and with non-canonical 
power - propagator for G^. This procedure can be immediately generalized to any number of 
loops. Furthermore, one should realize that the BT procedure, in short raising powers, does not 
introduce singularities through B{xg), a part from the singularities in the external parameters 
of the original diagram. Therefore, before performing the Xg-integration we move the integration 
contour into the complex hyper-plane, thus avoiding the crossing of apparent singularities, see also 
the work of ref. 0- 

As we have anticipated, in this paper we extend the minimal BT-approach to cover all two- 
loop, two-points functions. In dealing with the sunset topology we have discovered a remarkable 
property that will be generalized in this paper to all two-loop two-point functions but that remains 
true for all two-loop three-point planar topologies. For all these diagrams, essentially, we do not 
have enough external momenta to make the matrix H of Eq.(^ non-singular. As a consequence, 
a change of variables is always possible in the quadratic form V that makes the 5-coefficient of 
Eq.(|^ independent on Xg, at least if we use only one iteration of the algorithm. Therefore, in 
these cases, no distortion of the integration path is needed. In this way we are able to be closer 
to the original idea of the BT-method. However, as already stressed in the original paper [Q, 
this B will vanish at some threshold. We are able to show that this occurs at some non-leading 
Landau singularity of the diagram where additional analytical work is needed before starting the 
numerical evaluation. 

At the same time we have considered the problem of an efficient evaluation of infrared divergent 
multi-loop diagrams. Although the massless world is efficiently treated in QED and QCD within 
the analytical approach, any multi-loop calculation in the standard model or beyond is plagued by 
infrared divergences and we must contrive a scheme for dealing with them in a purely numerical 
way. For some of the two-loop two-point functions, the infrared divergent conhgurations are simple 
enough to allow for BT treatment with a consequent and straightforward analytical evaluation. 
The remaining topologies, in particular the one containing overlapping infrared conhgurations, 
require a novel approach. 

flTfi 


Our solution is derived by adapting the general algorithm of ref. 


the residues of the 


infrared poles and the hnite part of a multi-loop diagram can be cast into a form which allows 
for a reliable numerical integration. The algorithm is perfectly dehned for all integrands that may 
vanish only at the boundaries of the parametric space. Here the strategy differs somehow from 
the BT approach. In the latter the main attempt is toward raising powers in the integrand while 
here the powers remain untouched and the whole idea is about factorizing the singular behavior 
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into simple factors that can be integrated analytically, leaving non-singular terms to be treated 
numerically. 

The outline of the paper will be as follows: in Section § we introduce the dehnitions and the 
basic properties for arbitrary two-loop diagrams. In Subsection ^ we discuss the so-called re¬ 
ducible two-loop diagrams while, in Section ^ we present special cases of the general BT functional 
relation. In Section ^ we start the evaluation of two-loop two-point functions with four internal 
lines. In Section ^ a general class of two-loop diagrams, those that contain a self-energy insertion, 
is discussed. Explicit methods for evaluating diagrams with four internal lines are given in Subsec- 
while their Landau singularities are computed in Subsection 5.3. Further rehnements 


tion 5.1 


for diagrams with four internal lines are discussed in Subsection |5^ - 15.8| . The on-shell derivative 


of these diagrams and their infrared poles are analyzed in Subsection In Subsection |5.10| the 
tensor integrals are studied. The two-loop two-point diagrams with Eve internal lines are presented 
and discussed in Section |^. In particular we discuss in Subsection o the simplest case in this 
class, in Subsection its Landau singularities and in Subsection its explicit evaluation. The 
corresponding on-shell derivative and its infrared poles are shown in Subsection M- The most 
difficult topology with Eve internal lines is presented in Subsection [f.5|, its Landau singularities in 
Subsection 0 . its evaluation in Subsection 0-0 The corresponding on-shell derivative and 
its infrared poles are discussed in Subsection |7.9| while in Subsection [7.10| we analyze the tensor 
integrals of this class with a speciEc and realistic example shown in Subsection [f.ll . Finally, 
numerical results and comparisons for those few cases where analytical work has been done are 
shown and discussed in Section ||. Technical details are discussed in various Appendices. 


2 Arbitrary two-loop diagrams 

In this section we present basic dehnitions and properties of diagrams. An arbitrary two-loop 
diagram has the following expression: 

(o:, Tn\^ • • • : hi 1 0^7 ^^CK+l; • • • ; 75 )ma 7 hl 2 |/^7 7 • • • : ^^a-t- 7+/3 7 hs) 

2 e ^ a a+'y a+'y+P 

— (Fqi(rq 2 \{{k‘l + m^)-^ H (^1 + "^?)“^ IT + (5) 

^ i=l j=a+l Z=Q!+ 7 -|-l 

where u = 4 — e, with n being the space-time dimension, and where a, (3 and 7 give the number 
of lines in the qi, q 2 and qi — q 2 loops respectively. Furthermore we have 

ki = qi + E7=i VijPj, z = 1,..., a 

ki = qi-q2 + E7=i VijPj^ i = a -Fl,..., a -F 7 

ki = q2 + E7=i PijPji i = a 7 -h 1,..., a + 7 + / 3 , 

N being the number of external legs, 7 “ = ±1, or 0 and {p} the set of external momenta. 
Furthermore, /i is the arbitrary unit of mass. Diagrams that can always be reduced to combinations 
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of other diagrams with less internal lines will never receive a particular name. Otherwise we will 
denote a two-loop diagram with 

( 6 ) 


where G = S,V, B etc. stands for two-, three-, four-point etc, and a, (3 and 7 give the number of 
lines in the gi, q 2 and qi — q 2 loops respectively. Family of diagrams with the same number N of 
internal lines will be denoted collectively by S^, etc. 

Next we recall few basic properties of two-loop diagrams |^. In any two-loop diagram there 
are three one-loop sub-diagrams, called ay, /dy and a/? sub-diagrams respectively. 


Definition 1: the a/Jy diagram is overall ultraviolet (UV)convergent if a - 1 - - 1 - y > 4, logarith¬ 

mically divergent if a - 1 - y = 4, linearly etc. divergent if a -|- /? - 1 - y = 3,.... 


Definition 2: the a(3 sub-diagram is convergent if a-\-(3 > 2, logarithmically divergent if a-f/3 = 2, 
linearly etc. divergent if a -|- /3 = |,.... 

If any of the one-loop sub-diagrams diverge we have counter-terms associated with them. 
Therefore, in addition to the a/Jy diagram we will also consider the subtraction diagrams of 

Fig. m. 

etc., (7) 

The arbitrary diagram has the following representation in parametric space: 


^Q /37 _ _ 




TT 


F(a + /9-l-y-4-fe) / d[x\ / d[y\ X{l-X) 


2 a 7 e/2 Q,_|_.y_3^£/2 

Vi 


X 


Q!+/3-|-7 a+q+'y 

2/1 + X! (-^2« + '"^1) Ul-OL-y+l — {Px Vl + X! Vl-a-y+lY'^ 

/=Q:+7+l Z=Q:+7+l 


4—a—/3— 7 —£ 


. ( 8 ) 


where the integration measures are 


1 a +7 a+y i /3+1 /3+1 

d[x]= n d[y]= Yl dyi6{J2 
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i=l i=l ® i=l i=l 

and parameter dependent masses and momenta are defined by 

R^-K^ „ Pi Pi; 




.Y(l - .Y) 


P = 

) X 


12 


1 -X X ’ 


a+y 


K — Pi + P12, Pi — ^ Xi Pii, P12 — ^ Xi Pi2i, 

i=l i=(y.-\-l 

a Q:+7 

= ^Xi {P^ + (Pui + 

i=l i=OL-\-l 

N a+7 

Pai = J2vtkPk, x= Y. 

k=l 1=0+1 


(9) 


( 10 ) 
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diagrams. 



etc. 


of Eq.(D and one of the associated subtraction sub- 


2.1 Reducible two-loop diagrams 

There are two-loop diagrams that can be easily reduced to products of one-loop diagrams, 
some of them occurring in the reduction of tensor integrals. A typical example is represented by 




dP-qi (Pq2 


){ql + ml){{q2+py + ml) 


Ao(mi) Bq , 


where Aq{Bq) is the one-loop, one-(two-)point scalar function |jT^. This function satisfies 


( 11 ) 


AB{p^] mi, m2, m3) = AB{p^-, mi, m3, m2) 


( 12 ) 


3 Minimal BT algorithm 


For an arbitrary diagram the application of minimal BT algorithm requires distortion of the 
integration hyper-contour. First of all let us discuss how this can be done: we rely on the fact 
that applying BT functional relations to does not introduce new singularities. Consider any 
diagram: the integration hyper-contour can be distorted away from its original real location, a 
simplex described by Xj > 0 and < 1, and when the possibility of this distortion ceases, 
we encounter a singularity of the function The difficulty in locating the singularities lies in 


imagining what happens in the multi-dimensional (complex) space of integration. Consider now 
a simple example 


G = 


[ dx f dyV ^ "{x,y)= f dyf{y), 
Jo Jo Jo 


(13) 


V{x,y) = H x"^ + 2 K{y) X + L{y). 
Minimal BT algorithm transforms G into 


[1 + ^lnV{x,y) 


rl rx 

G= dx dy . 

Jo Jo o[y) L 

^ ^ dyTr^{[^ + X{y)]\nV{l,y)-[y + X{y)]\nV{y,y)}, 
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(14) 


(15) 
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with b = L — K"^/H and X = K/H. Let hq be a root of 6, i.e. &(|/o) = 0. What happens to /(i/o)? 
Since &(?/o) = 0 we can use another identity 


1 + 


2 (1 + e) 


{x + X)d^ V ^(x,yo) = 0, 


(16) 


to derive 

f(yo) = - r,n\ , [ dx{x + X)d^V{x,yo). (17) 

Z\L + e) Jyo 

If there is a pinch singularity a.t x = Xq with yQ < xq < 1 then {xq, yo} is already a singularity of 
G. If not we integrate by parts and obtain 


(1 + 2 e) f{yo) = {yo + X) V-^-^{yo, yo) - (1 + X) l"-i-^(l, yo). (18) 


If V{yo, yo) = 0 or 1/(1, y^) = 0 then x = 1, ?/ = ?/o or x = ?/ = ?/o is an end-point singularity of G, 
if not then no new singularity is introduced since f{y) is analytical in y = yo. 

To understand when minimal BT eventually fails we consider another simple example which, 
however, is general enough to complete our discussion. We start with a quadratic form 

V{x,y) = Hx"^+ 2K{y)x +L{y), (19) 


H = ml 


K = ^{s-ml-ml)y L = mly'^ + {M^ - s) y + (20) 


This is actually a Gq function corresponding to a vertex correction for W (off-shell) — ud. Applying 
minimal BT gives a b{y) which is a quadratic form in y. The two roots b{y±) = 0 are coincident 
when 


s = 


2 ml 


mlhf; - MlM'i + ± M", mj)!, (21) 


which is the well-known anomalous threshold. Here A is the usual Kallen-function |H|. Let y^^^ 
the solutions of H = 0, which are now branch points in Eq.(^). These roots are coincident for 

X = xth = (ml Ml-2mlMl-mls + mlMl-mls-Mls + s^) A"^(s, ml, ml), (22) 

and one can easily verify that with x and s given by Eq.(^^ and Eq. (0 respectively we have 

y± = yL = yR. (23) 


In this case, if Xth G [0,1] and y± G [0,1] we cannot distort anymore the y integration contour 
to avoid y = y± without crossing a cut of the logarithm. However we can still go on with our 
procedure if we distort the x integration contour in order to avoid x = Xth- This will be possible 
unless, for some specihc conhguration of the external parameters, we have Xth = 0 or xth = 1- 
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In our particular example this would require an unrealistic relation between W, Z and up, down 
masses. Even in this case we still have a chance of continuing the derivation. Indeed we have used 
logarithms of quadratic forms which give a pair of branch points which can pinch the integration 
contour exactly where we would like to deform it. However, when only logarithms are present and 
not logarithms squared, we split the logarithm as 

lnml + rj(y-y^,y-yj,)+ Hy-Pi), (24) 

i=L^R 

where y is the usual Veltman’s function, p(a, h) = ln(a6) — In a — In 6, and logarithms are dehned 
with a cut along the negative real axis. The hnal expression for the diagram will be the sum 
of terms proportional to r^-functions and terms of the form ln(|/ — yi)/b{y). The former do not 
create a problem in distorting, when needed, the integration contour. For the latter, however, if 
one of the roots of b, or both, is real and internal to the interval [0,1] then the distortion for the 
corresponding term must be examined with care. If the imaginary part of yi is positive (negative) 
then we move the contour into the negative (positive) imaginary half-plane, so that the cut of the 
logarithm will never be crossed. If yi is real and negative no problem will arise and we can distort 
in any of the two ways. If its is real and positive then the distortion is hxed according to the i 6 
prescription. The reason why we cannot apply this argument to cases where we have the presence 
of In^-terms is simple: after the splitting of the logarithm we will encounter terms proportional 
to ln(?/ — y^) ln(j/ — j/^). When the roots are complex conjugated we will have cuts both in the 
positive and in the negative imaginary half-planes and the possibility of distorting the contour 
ceases when they approach the real axis. Note that these quadratic terms are always present for 
ultraviolet divergent diagrams. 


4 Special BT relations 


In this section we present some explicit examples of the Bernstein-Tkachov functional relations. 
Let V{x) be a quadratic form of a vector x, 


V{x) = x^ H X + 2 X + L, 


(25) 


then the solution to the problem of determining the polynomial V is as given in Eq.(^). There are 
special BT functional relations: typically, when we have a quadratic in two variables y, z of the 
form V = hz^ + 2 kiz -|- 2 k 2 y -l- I then the matrix H is singular. Here we use 


1 - 


y 


d„ - 


y \ 2 (/i -1- 1) 


\ h) ^ h 


For a quadratic in one variable, V = hy"^ + 2 ky + I we use 


2(p + l) 


^ k\ d 
^ ^ h) dy 


pA/2+l ^ 


hl-k'^ 

h 




(26) 


(27) 
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while for 1 / = ky + I we use 


1 d ■ 
/i + 1 ^ dy- 






1 1 - 1 / di 

y + 1 k dy- 


^/^+i 


V^. 


(28) 


5 The ^4 family 

There are two diagrams in the S' 4 -family, both with a = 7 = 1 and with (3 = 2. The hrst is 
and is given in Fig. It is overall logarithmically divergent and the ay sub-diagram is also 
logarithmically divergent. We have another diagram with the same number of internal lines, given 
in Fig. This second diagram is given by the following integral: 



Figure 2: The two-loop diagram of Eq.(^). Arrows indicate the momentum flow. 




dJ^qidJ^q2 


{ql + ml) ( (gi - + "li) (gi -F ml) {ql + ml) 

Therefore, this case is trivially reduced to known diagrams. If m 3 7 ^ m^ then 

1 


(I;mi|l;m2|2,m3,m4) = 


2 2 

ml — m 3 


i clll ,,2e 
7Tb = fi 

otherwise we will have 


mi, m2, m3) - mi, m2, m4) 

d^qi d^q2 


{ql -F ml) ( (gi - g 2 + p)^ + mfj {ql + m|) 


(I;mi|l;m 2 | 2 ,m 3 ,m 4 ) = (e - 1 ) S' 33 ( 0 ; mi, m 2 , m 3 = m 4 ), 
where 5*33 has been introduced in ref. as the following integral: 


vr^ S'33 = 




2 e 


d^qi d^q2 


e — 


1 J (?i + mf) [{<]i - <]2 + pf + mfj 


(29) 


(30) 

(31) 

(32) 


(33) 


















Figure 3: The second diagram of the S 4 -family, Eq.(|^, which is evaluated in terms of (ref. |^) or 
533 diagrams (Eq.(|3^) at zero external momentum. 


The topology is given in terms of the following integral: 


4 0121 _ ,, 
TT O — 


2e 


d^qi dP'q2 


( 9 i + ”>1) ((91 - 'nf + >"2) (92 + ™i) ((92 + pf + ini 


(34) 


One way to evaluate starts with the introduction of a Feynman parameter x for the qi- 
loop; successively, we integrate over qi and, finally, combine the g 2 -dependent propagators with 
parameters y, z. After performing the g 2 -iiitegration we obtain 




X 


^5, 


{pi - 


r (e) / dx dy 


dz 


a;(l - a;)] (1 - y) 


se/2 


X 


-l-e 

121 


+ 


2 -e/ 2 ,_, 


X 


121 


where the quadratic form Vi is dehned by 

X 121 Vi z) =z^ - /i ^4 z + {fil- yl) y + yl, 
and where we have introduced additional auxiliary quantities. 


(35) 

(36) 

x(i-x) x(i-xy ™ ‘ 

In Eq.f^Tl) we have assumed that = —s and s > 0. Scaled masses are, therefore, defined to be 
= y? s. If, instead, is positive, we will have to change the procedure as described in the 
next subsection. 


2 _ (1 - x)yl + xyl 


/^?2 


2 _ 1 I 2 _ 2 

^■34 ~ + A'-S ^^A■ 


5.0.1 The case = —t with t < 0 
Starting from the quadratic form 

t/ = (y — z) + z + (1 — y) + z{l — z) p^, (38) 

with p^ = —t and t negative we redehne scaled masses to be t and obtain U = — f x, 

X = + yl^ z + {yl- yl) y + yl, with yl^ = l + yl- yl. (39) 

The signs in this quadratic form will have an effect on the raising procedure, as described below 
in Subsection 
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5.1 Evaluation of method I 


In the following sections we will be able to present several derivations for and to compare 
the corresponding nnmerical resnlts. In order to evaluate we may use the following relation: 


X 


-l-e 

121 


4 r 1 1 - 

{x, y,z) = - -— 1 + - {Pydy + y, z), 


121 


with polynomials Py^z given by 


Py Pz 2 2 ^121 ^34 ^ ha 

If = —t and t is negative, the dehnition of will change accordingly, 


^121 


h 43 + 4/ia; > 0, 


(40) 


(41) 


(42) 


so that 6^21 is never zero in the t-channel, as expected. 

We will also use an additional quadratic form, X^ 2 i) which is related to x^ 2 i (^0- 


),by 


^121 h, z) = x(l - x) X ^21 (^, y, ^) 

= x{l - x) -z{l - z) + pIz + pl{y - z) + {1 - x)pI + xpl {1 - y), (43) 

The next ingredient needed in the evaluation of the diagram is integration by parts which is 
performed with the help of the following relations: 

dz z'^dz In"" X^2i y^ = y'^ -^121 v) ~ ^o,n In™ X^2i 0 ) 

dzz^-^ In™A,,,, 

/o 


/o 


n 


dz z" ^ In™ A .21 {x, y, z), 


(44) 


^ dyy"(l In” Ay, (i, !/, j) = -j”(l - z, z)) 

- j' dy[ny"-\l-yrl^--"-y’\l-yyl^-'] ln”Ay,(i,s,y. 

Furthermore, we will use 

%»”(! - y)-'^-' \xrX^^y{x,y,z) = -{l- zf/'^ ^"*^,,( 1 , l,z) 

J z 6 

+ / <ij/(l - A , (i.s.z) - lii"*A (i, l,z) 


(45) 


(46) 


Before discussing the hnite part of we consider its ultraviolet poles. As usual, the pole parts 
can be computed analytically. The double pole of is easily derived and the corresponding 
expression, where in the spirit of the MS scheme we include a bunch of constants, reads as follows: 

2 


0121 

double—pole 


2 uv^ 


(47) 


10 














where, for the constants, we obtain 


A 


uv 


7 + hiTT + In 4 t- 


(48) 


Here, 7 = 0.577216 is the Euler’s constant and /i is the arbitrary unit of mass. The single-pole 
requires some extra work to be discussed in Subsection pA| . 


5.2 Single-pole and finite part of 

For computing the single ultraviolet pole of we use the relation 


hx - 


X 


= 1 - {zd^ + ydy) lnxi2i> 


121 


and derive the following result: 


Sl^iie-poie = ^ + constants. 


= 2 ( 7 -f luTT-f In - 1 + 2 [ dx [ dy \nx,r,Jx,l,y). 

Jo Jo 


where we also include a bunch of constant into the dehnition and where 

nl 

^;)-l + 2/ 

Since X121 (^. ^,y) = - Av + y-h *'* find 


(49) 


(50) 


(51) 


(52) 


which is the hnite part of the one-loop two-point function corresponding to internal masses m 3 
and m 4 . This is exactly the result expected, from general arguments, to compensate the over¬ 
lapping divergency coming from the associated subtraction diagram of Fig. This result proves 
cancellation of poles with logarithmic residues, leaving only harmless poles, i.e. those with as 
residue a polynomial of hnite order in the external momenta. Summarizing our Endings for the 



Figure 4: The subtraction diagram, containing a one-loop counter-term (represented by a x) associated 
with the two-loop diagram of Fig. 
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(53) 


ultraviolet divergent part of we obtain 

>5siniie-poie = ^ (2 ^uv “ 1) + + 2 (^ - ^ dy InXi 2 i (x, 1, y), 

where one should note that Xi 2 i (^5 x-independent. 

For the hnite part of we have to consider the following quadratic forms that arise after 
integration by parts: 

Xi 2 i = {y- + zyl + (i - y)yl - z{i - z), 

Xi 2 i y^ y) = yyl + (i - y)y'l - yi^ - y)^ 

Xi2i y^ 0) = yy^l + (i - y)yi^ 

Xi 2 i 1) 2 /) = (1 - 2/)h3 + yyl - 2/(1 - y)- (54) 

Before continuing the evaluation of the hnite part of 5^^^, however, we need to understand its 

analytical structure. The tools for extracting all the singularities of the diagrams are, as usual, 
represented by the corresponding set of Landau equations. 


5.3 Landau equations for 

Before starting the evaluation of this topology it is important to know as much as possible 
about its singularities as a function of and of the internal masses. The corresponding Landau 
equations are given by 

ai{ql + ml) = Q ^2 ((gi - ^ 2 )^ + = 0 , 

«3 (gi + m^) = 0 0(4 {{ql + pf + ml) = 0, 

and also by 


o^igiM + '^ 2 (gi ~ — 0, —Q: 2 (gi ~ + c^^(l2^J, + ct4 (g 2 + p)/i — o. (55) 

The leading Landau singularity occurs for Oj ^ 0,Vi We multiply the two equations Eq. (|55|) by 
ql^l, q 2 fi and respectively. This gives an homogeneous system of six equations. If all are 
different from zero, the singularity will occur for 


g? = 


-mr 


ms 


qi-q 2 = \ (ml - ml 
The Landau equations become as follows: 

—2 ml cxi + (m3 — ml — m\) a2 = 0, 


g 2 ^ = 


-m 


3’ 


P ■ q2 = \ {s + ml 


ml). 


/ 2|2 2 \ i /2 2 | 2 \ n 

— (m3 + m^ — m2) «i + (m3 — m^ + m2) 0:2 = U, 


P ■ gi «i 

.2 I .^2 


p ■ gi + - (m4 


mg 


s) 0(2 = 0 , 


(ml + ml — ml) 02 + (m^ — ml — ml) 03 + 2 p ■ gi + ml 


m. 


mn 


0(4 = 0 , 


m 


mn 


m3) 0(2 — 2 m3 Os + (s 


m. 


m. 


0(4 = 0 , 


s + ml — ml — 2 p ■ qi 02 + (s + mg — ml) 03 — (s — mg + ml) 04 = 0. 


(56) 
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There are two compatibility conditions for the first three equations that can be used to derive 


p-qi = 


mi 


2 (mi + m 2 ) L 

Furthermore we must have a relation among masses, i.e. 

m\ = (mi + m 2 )^ 

From these considerations we derive that 


s + (mi + m 2 f - m] 


m 2 

tti — -0:2, 0:2 — -7-;-r 

mi m2 (mi + m2) 


|-(mi + m2)^ 0:3 + ^ s - (mi + m2)^ - m4 q;4|, 


is a solution for the system if the following condition is satished: 

s = (mi + m 2 ± m 4 )^. 

Note that we have obtained the leading Landau singularity, i.e. all 7 ^ 0, as 

s = (mi + m 2 ± m 4 )^, m 3 = (mi + m 2 )^. 


(57) 

(58) 

(59) 

(60) 

(61) 


Naturally we have that 

s = (mi + m 2 ± m 4 )^, 0:3 = 0 , (62) 

are sub-leading Landau singularities corresponding to the normal and pseudo thresholds of the 
reduced diagram where the line corresponding to the propagator -|- m| is contracted to a point. 
Similarly, 

s = (m 3 ± m/iY, oi = 02 = 0, (63) 

are the sub-leading ones corresponding to the reduced diagram where the qi sub-diagram is shrunk 
to a point. 


5.4 One iteration for 

Having discussed the analytical properties of the diagram we continue the derivation of an 
integral representation which is particularly suited for numerical treatment. After one iteration of 
the raising procedure the result for is still very compact. First we introduce a combination 
of logarithms. 


L±{x^ y^z) = ^ lnX(a;, ?/, z) ± |^ln(l — y) -|- In a; -|- ln(l — x) 
and obtain for the hnite part of 


(64) 


Cl21 _ 
^fin — 


H 121 ^0 


dx 




dy J dz + 


/o 


dyip^+n^^ 


/■I /■! fy r 1] 21 1 

^ Jo Jo Jo ~ T ~ 2 
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( 66 ) 


where B follows from Eq. (|4l|) and becomes 


^121 = ^(1 “ ^21 = ^(1 “ /^34 - 4/i?2, 

with /if 2 /i |4 dehned in Eg. (p7D . and the varions integrands appearing inside Eq.(|65[) are as 

follows: 

/ 3 I 21 = - In X 324 (x, //, z) {4 (/if 2 -Xz^)-2 L+{x, //, z) - Xz {b z - } 

dnX^^i {x,y,z) 


-A{Xz‘^- /ifa) (- 


2 /- 1 


4 


(67) 


J121 = _4 (X//2 -/if^) \nX^^^{x, l,//)L_(x, 1,//) - (p ^4 + 2//) (X//^ - ^^ 3 ) 

X \nX^^^{x,y,y)L+{x,y,y)- \nX^^^{x,y,^)L+{x,y,^), 


( 68 ) 


1 2 ^ 
/(•a = - In X (- X - 2 3 + 55 -’f. 


(69) 


Here (C(z) is the Riemann zeta-fnnction [|Tg, X = x{l — x) and two additional qnantities have 
been introdnced: 

^■12 ~ (1 ~ hi + ^h2) h34 ~ 1 4" hs ~ h4- ('^O) 

The ‘+’-distribntion is dehned, as usual, by its action on a generic test function f{x)\ 

(n) 

If we stop the chain of iterations at this level we will have the advantage of a result that contains 


only one denominator, namely 6321 - ^ second iteration, described in the Section makes the 
total integrand smoother at the price of introducing additional denominators. If = —t with t 
negative we dehne scaled masses as = —fx^t and derive 

J 3 I 21 = {-4 {Xz^ + yl^)+ 2\x{bz^ - yl^z) L+(x,//,z)} \nX^^^{x,y, z) 

, o 9 X (x, y, z) \ 

+ 4 (Xz^ + yl^ ’ 

= {Xy^ + /ifa) (hL + 2 y) In X 324 (x, y, y) L+(x, y, y) - yl^ yj^ In X 324 (x, y, 0 ) L+(x, y, 0 ) 


+ 4(Xy +yu) lnX 324 (x, l,y)L_(x, l,y), 

1 2’^ 

/I2 i=(_x + 2 ,,y 


(72) 


instead of Eqs. (§3)-(§^. As usual, the iterative procedure of raising powers in the integrand 
produces apparent singularities that make the hnal result unstable in known regions of the x- 
integration. In the next subsection we will explain how to cure these instabilities. 
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5.5 Distorting the integration contour for 


Obviously, there are numerical instabilities when the zeros of b. 


Liic Z.C 1 WD -ji <y^ 2 i (Eq.(P^), are real and 
internal to the integration domain. This can be avoided by distorting the x integration contour 
in Eq.(^) into the complex plane. However, this possibility ceases when the roots of the 
quadratic form (Eq.(|43|)) pinch the real x-axis. Given that X^ 2 i = ax^ + 6x + c, 

this situation will occur when we simultaneously have 


x° = x*! =-and 

+ 2a 


r = 4 ac. 


The condition for coincidence, i.e. x° = x° G [0,1], is 

/i34 = 2 (/Xi + ^ 2 )- 


(73) 


(74) 


The remaining two conditions, namely Imx° ^ = 0 and x° = x° = x° = x!J_, require 


(A-A)y + (~ -= 


(75) 


where = pi + /r 2 - There are, therefore, two solutions for given by 


z = 


h+ ± (h 



1/2 


(76) 


and, as long as we can distort the z integration contour to avoid these points then the x contour 
will not be pinched. This is always possible unless a new pinch will occur, which indeed is the 
case if /is = /i+. Note that // = 0 and z = /i+ is outside the physical region where 0 < z < //. 
By inserting /is = /i+ inside the relation /is 4 = 2 (/ii + ^ 2 ) we obtain /i+ ± /i 4 = 1 and the two 
correspond to the leading Landau singularity. 

There are also logarithms of quadratic forms in two variables, x and y. For (^5 v) 
x-contour is pinched for 


1/ = h+ + ^ {h+ 




(h3 - h+) (h3 


h+ 




(77) 


Therefore, the y contour will be pinched if /ig = /i^ or /ig = /i^ + 4 /i^.. Only the former is physical, 
corresponding to // = /i+, while the latter corresponds to y = — /i+. 

For I/) 0) the x-pinching occurs for y = —/i^/(/i| — /i^) so that it is enough to distort 

the //-integration contour avoiding this point. Finally, for 1,//) the solution is 


y = y+±{yl- (78) 

showing, again a pinch for /ig = /i+. Around these points we must introduce an alternative 
derivation. 
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5.6 Around the leading Landau singularity: method II 

When we are in the regions of external parameters that correspond to the leading Landau 
singularity of an alternative approach must be used. We easily derive the following relation: 

S'^^^ = (e — 1 ) J dx S 33 (^x‘^p‘^;mi,m 2 , , with =—p^x'^ + {p^+ ml — ml)x + ml, (79) 

where S'33 (Eq.(^)) has been introduced in ref. [0. The finite part of S'33 is 
(e - 1 ) S'fg" (xV; 1^1,1712, M 3 .) = 


- f dy f dz [ln^(|/,z) + 

Jo Jo L y z — 1 

M^ / M^\ 7 1 

‘■'V - 2 - 2 « 2 )' 


(80) 


(81) 


and ^ is given by 

^{y, z] /ii, p 2 , Ps) = zy{y - z) [(1 - z)x‘^ - /ig] + (1 - ^) [{pi - pl)y + pl 
and must be evaluated at /i| = Mj/ s. 

5.7 Two iterations for 

In this section we return to the method described in Section E3 an additional iteration of the 
procedure requires the consideration of new quadratic forms in two or one variables. In particular, 
when we apply integration by parts, some care is needed in dealing with terms containing ln(l —?/). 
For instance, we will use 

r-l ry rl rl 

Jo Jo ln(l - 2 /) = dzj^ dy In y, z)\n{l - y), 


10 


dz ( dyy^ In V, z) \n{l - y) = [ dz [ dyy^ In ^121 

Jz Jz ^iV) 


+ jo dz\nx^ 2 iA^^^A dy{l -y)'^\ny. (82) 


With only two variables we use 


fl pI 

dyy''lnx 32 iA^y^y)M^-y) = dyy'^ln- 


Xi2iA^y^y) 

X 121 


ln(l - y) 


1 1 
7 ^ ^ 7 ") 

n + 1 ] 

Xi2i (a^i2/i0) 


/ dyy"^ In {x, 2 /, 0 ) ln(l - y) = f dyy'^ In ln(l - y) 

JO jO ^121 


1 1 
r—T lnx,2,(x,l,0)^ -. 
-L j = l J 


(83) 
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The second iteration of the rasing procedure is based on the following set of identities: 

4 fl r n 

InXi2iI/) = -7 - [7; {y + yi)dyXi2ii^^y^y) l-lnXi2i(a^>l/>l/) 

+ X 121 y^ y) In X 121 y,y)}, 

8 1 

In^X121 I/) = 7 - {l^iy + ^1) dy[xi2i{^^y^y) -Xi2i(a^>l/>l/)lnXi2i(a^>l/>l/) 

1 1 

+ 2^121 y ^ y ) T121 I/)] - 2^121 (^’T121 ?/)}• ( 84 ) 

Furthermore, we will use an additional relation: 

/ dyy^dyX^^i (^’ T 121 y) = 1) In™ X 121 1) 

J 0 

- <5n, 0 X 121 (^’ 0’ 0) X 121 0> 0) 

-n dyy^~^Xui(^^y^y)^'^'^Xi2i(^^y^y)- (85) 

J 0 

Here we have introduced the following quantities: 

Yi = ~(l-yl +yl) , \{l,yl,yl), ( 86 ) 

with special cases given by X^ 2 i(^’ ~ /wland Xj^2i(^)0)0) = ^‘x- Similarly, for X^ 2 i(^’ 1)1/)) we 

obtain 

4 (1 

lnXi2i(a^) 1)//) = -7 -|^(l/ + l2)'9yX^2^(a;, 1,1/) 

®2;121 ^ 

X [l-lnXi2i(^’l’^)] +Xi2i(a^) l)l/)lnXi2i(^’l’^)}’ 

8 1 

in^ X 121 (a^) 1) y) = 7— {^ (y + 12 ) <92/ [xi 2 i (^) 1) //) - X 121 (a^) 1) y) in Xi 2 i {x, i, y) 

" 2;121 ^ 

+ ^Xi2i(a^)l)//)ln^Xi2i(a^)l)//)] - ^Xi2i(a^)l)l/)ln^Xi2i(a^)l)l/)}- (87) 

Furthermore, the following identities hold: 

/ dyy^dyX^2i (^’ l’^™ X 121 {x, 1, //) = x^2i (^’ 1> 1) 1'^'" X 121 (x, 1,1) 

J 0 

- 0 X 121 (^’ 1’ 0) l'^™ X 121 (a^) 1) 0) 

-n dyy^~\^^^{x,l,y)\n^Xui(^^'^^y)^ 

0 

/ dyy^dyX^2i (^’ 1’ = X^ (a^) 1) 1) “ '^n, 0 X 121 (^’ 1’ 0) 

J 0 

-nj^ dyy^~^Xi2i(^^^^y)^ (88) 
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where we have introduced a new F-factor, 

¥2 = -^ (1 - h4 + hs) > ^2;121 = /^4), (89) 

with special cases given by X 121 ~ hi) X 121 (^) 9) 9) = hi also x^ 2 i (^) 1) 9) = p|. Finally 

we have to consider terms containing X;^2i(^’h)9) which give 


In X121 h, 9 ) = X121 (^’ h, 9 ) In x^2i h, 9 ) - - X121 (^’ h, 9 ) 

rXi2i(a^)h,9) 


d, 


Xi2i(a^)h,9) 


3;121 


1 - X 121 (a^)h,9) 


a, 


, lnXi2i(^’h,9) 

"3;121 


ln^Xi2i(a^)h,9) = -2{--x^2i(^’h,9)ln^Xi2i(^’h,9) - 

rXi2i(a^)h,9) 


I-X 121 ix,y,o) 


a, 


2 L 


1 -Xi 2 i (a^)h,9) 


a, 


1 - Xi2i(a^)h,9) 

, lnXi2i(^’h,9) 

^3;121 

r^121 (^’ h) 9 ) 2 / p,\i t 

In Xi2i(a^)h,9)J|, 


a 


X 121 (a^)h,9) 


3;121 


3;121 


where the special case is X^ 2 i(^’h) 9) = yfi^ + (1 — l/)hl- For this combination we use 


/ d|/|/”<9yX^2i(^’h,9) = Xi 2 i(^’1’9) - (5„,oXi2i(a^)9,0)-n f dyy^ 
a 0 a 0 


Here we have dehned 
For terms involving X 121 (^) h) the second step will be as follows: 


h _ 2 _ 2 

^3;121 h3 ha;- 


inxi2i(a^)h,^) = -7 — {ydyXi2iix,y,z)[^-'^^Xi2iix,y,z) 


'121 


1 


+ -(^ + >2) <9^X121 h) [1 - In X121 h, z) 

+ Xi2i(a^)h,^) lnXi2i(^’h,^)}, 


(99) 


dyy^dyX^2i (^’ h, 9) In”^ X 121 (^’ h, 9) = X 121 (^’ 1’ 9) In™ x^2i (^’ 9) 

- 5n, 0 X 121 (^’ 9, 0) In™ x^ 2 i (^) 9, 0) 

-n f d|/|/”"^X^2i(^’h,9)ln™x^2i(^’h,9), (91) 

a 0 


(92) 

(93) 


18 



2 1 
Xi2i = I —{ydyXui (x, y, z) [i - in Xi2i {x, y,z) + - in^ Xi2i y^ 


121 


+ ^ (^ + >2) <9^X121 y,z)\\--\n X121 {x:y^z) + ]^ in^ X121 y^ 


Integration by parts can now be performed according to 
ry 

/ dzz^d,Xi 2 i y^ X121 = y^ ^121 y) ^121 y) 

J 0 

- <5n, 0 X 121 (^’ 0) ^121 (^’ 0) 


(94) 


n / Xi 2 i(a^>d,^)ln'"Xi 2 i(a^>d,^), (95) 


p fy _i 

/ d 2 : 2 ;'^( 9 ^Xi 2 i y,z) = y'^ X 121 d, d) - <^n, 0 X 121 (a^, d, 0) - n / X 121 d, 2 :) • (96) 

0 ^0 


Furthermore, we exchange to order of integration, 

rl j-y rl rl 

/ dy dzdy = dz dydy, 

Jo Jo Jo Jz 

to perform the remaining integration by parts according to 

^ dyy^dyX^ 2 i (^’ d, z) In”^ x^ 2 i (^> d, z) = X 121 1 , z) In™ y^^i (^> 

- X121 (a^> ^) In'" X121 (a^> ^) 

- n ^ di/i/”" V 121 (a^> d, ^) In’" X 121 (a^> d, ^) 

^ dyy^dyX^ 2 i (^’ d, ^) = X121 (a^> 1 > ^) “ X121 (a^> ^) 

-n^ dyy^-^XuA^^y^^)- 


(97) 


(98) 


To summarize, the following quadratic forms will appear in the hnal result as arguments of the 
logarithms: 

77 i 2 i (X: y, z) = x{l - x) [{y - z)yl + zyl - z{l - z)\ + {1 - y) [(1 - x)yl + xyj 
77,21 {x, y, y) = x{l - x) \yyl - y{l - y)] + (1 - y) f(l - x)yl + xyj 


77^21 (^’ d, 0 ) = x{l - x)yyl + {I - y) [(1 - x)yl + xyl 
77^21 {x, 1 , y) = {l- y)y,l + y^l - y{l - y). 


(99) 
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Collecting the various terms we obtain a final answer that can be cast into the following form: 




*1: tf ^ s t^if 

3 “'0 “'0 n =0 ^0 n =0 '^0 ^0 n=l n=l 


( 100 ) 


Since we are not dealing with an analytical result and given that the number of terms in Eq.( | 100 |) 
is O ( 1.6 X 10 ^) we will not present its final explicit form in this paper. 


5. 7.1 Zeros of b-functions 

In the final result for we encounter negative powers of four different 6-functions whose zeros 
are as follows: 


- 6^21 of Eg. ( 1 ^). There are two roots for x, i.e. 


= TT 


h34 




( 101 ) 


where = 1 -1- p.3 — /i^. These roots are real and internal to the interval [ 0 , 1 ] if < 
4 (/xi - /ia)^ or /i^4 > 4 (/xi /is)^ 


64.421 of Eq.(^6[). There are three possibilities: 

1. /ii + — h4 ^ 1) therefore 64,424 can never be zero. 

2. (1 — /i4)^ < (/ii -|- ^2)'^ < (1 + h4)^i when there are two values of x where 64.424 = 0, 

1 


^ 1 + — 
x^ — 


2(1 + ^4) 


(1 + /i4)^ + hi ~ +2 ’ ^1’ ^^ 2 ) 


( 102 ) 


3 . (yUi -|- /i2)^ < (1 — h4)^) when we have four values of x where 64.424 = 0 . The new pair 
of points is given by 


= 


1 


(1 - - iii ±((1 - . A i*!) 


( 103 ) 


^ 2;121 


2 (l-h4)^ 

of Eq. (| 89 |) . This is x-independent and vanishes for s = (m3 — m4)^ or s = (m3 -|- 


The latter corresponds to the normal threshold relative to the two-particle cut of the 


diagram |T3 
b 

3;121 


6., of Eq. (P 3 D . There are two roots for x, i.e. 


x± = 


2 hi L' 




( 104 ) 


These roots are real and internal to the interval [ 0 , 1 ] if /i| < (/xi — H2Y or /i| > (/xi -|- /i2)^- 
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5 . 7.2 The logarithms 

To discuss the logarithms and their branch points we write 

lnX^ 2 i Vy = ln(aa;^ + bx + c — i S), (105) 

a = z(l - z) - (y - z)/U 3 - z/ul b = -a + (1 - y)(/U 2 - Mi), c=(l-y)/ul (106) 
Let us denote the two roots of the quadratic by There are three sub-cases to be discussed, 

- y, y) where a = y{l — y) — yyl and with roots 

- Xj^ 2 i 1 ) y), where a = y{l — y) — {1 — y)M 3 ~ yhl, ^ and c = 0, with roots x\ 

- X^ 2 i (^) y, 0 ) where a = —yyl, and with roots 


5.8 Evaluation of methods III and IV 

can be cast into different forms, all suitable for numerical integration. Let us consider 
again Eq.(|S3|) and change variable, y = 1 — y'. We obtain 


5^21 = _2 


TT S 


T (e) / dx dy 


nl-y 


dz 


x(l — x) 


1 - 6/2 


y /2 


X 


{mI - 


y ^ + 

^121 




"121 


where the quadratic form y^ 2 i is now given by 

Xi2i ^-y,z) = z'^ - Muz + {mI - MDy + mI- 

At this point we apply a special differential relation: 


y 


1+ - dy+ — {a 


(Mx - Ms) X 
1 

■ M 34 ) 


,,-1-e 


1 


121 


= --dyX 


—e 

y Ai_ 2 i’ 


- 1-6 


Y ^ = - b' Y 

A121 4 121 A121 


e “ 2 e ' 2 

where now b'_^^_^ has the advantage of being x-independent, namely we get 


^121 Ms, Ma)- 


(107) 


(108) 


(109) 


( 110 ) 


This second method is closer to the original BT approach in the sense that the denominator, 
after the raising procedure, is a function of the external parameters but does not depend on 
residual Feynman parameters. The zeros of b' , Eq.(|110|), are the normal and pseudo thresholds 


corresponding to the two-particle cut of Fig. ^ They are non-leading Landau singularities. The 
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fact that is x independent means that, with this method, we can perform the a:-integration 
without having to distort the integration contour. After integration by parts we obtain 




/o 


h' 

121 

1 .1 

dx I 
0 


1 r /-i / r^-y 

dx / dy 
Jo 


10 


dzJS' + Ji?)--4-- 

+ ^20^) + g + 4 }’ 


0 


( 111 ) 


where the various integrands are as follows: 


4i = 2 (/i 3 -^ ) lnx^ 2 i(^’l~ 2 /’^)“( 3 h 3 + h 34 ^- 5 ^ ) \n I - y, z) C-{x,y, z), 


7'121 _ 
"'21 — 


11 1 
“2 ^ 34^3 + 2 ^34 + hs “ 1 + (i- “ hL + hD 1/ + (1 + 2 ^34) ~ 

1 

2,.2 


X In X ,21 1 - 2/) 1 - y) I/, 1 - 2/) + o /^3h34 InX 121 1 “ 2/, 0) >^-( 0 ;, 1 /, 0), 




2r^or-o. - ..121' 

1 / lnxi2i(3:, l-y,zy 

2 I y 


Jlo^ = o 1 - 2/) 1 - 2/) 1-2/)+ lnXi2i (^>1>2/) 1,2/) 


2 L 


( 112 ) 


In Eg .( 11121) we have introduced special combinations of logarithms. 


^^-{x,y,z) = ln?/-lna:-ln(l-a:) - y, z), 

C+{x, 1 , 1 /) = lnx^ 2 i(^>l> 2 /) + lna: + ln(l - a:) -ln(l -//). (113) 


If requested, the procedure can be iterated without introducing numerical instabilities. However, 
the method fails at the non-leading Landau singularities corresponding to the one-loop sub-diagram 
with the largest number of internal lines. Therefore, in these regions we have to modify the 
procedure. Here, = ±2 /is and x^ 2 i = {z ^ ~ hi) 2 /j correspondingly we change 

variable, z' = z ^ and obtain 


lo Lk ±^*3)^) X'i 2 i ‘ + X121 = - h)!/■ (114) 


2-V2,_. 


After a new change of variable, z"^ = t, we use M{t, y) = t + (/x^ _ ^ 2 ^ y (derive 

{(;■» - l4) ; 1 } M-'-‘ = ; a.,} ^biM - I In" M + O [( 


(115) 


and arrive at the hnal result. For the diagram we have also another method (IV) which makes 
use of the following general property: 





(116) 
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Figure 5: The two-particle cut of diagram of Fig. 
After integration by parts we obtain: 


= - {^ dx dy dz 


I/- 1 


d i 'I 

~ Jo Jo +2 + 


3 1 


with Li given by 


Li{x,y,:) = lii(l - a) - hi a: - lii(l - x) - hix,,, {x,y,x)- 


(117) 


(118) 


Method IV represents the simplest results for but it is not generalizable to more complex 
topologies. 


5.9 Derivative of and infrared poles. 

The approach described in this paper is primarily intended for evaluation of massive multi-loop 
diagrams. However, QED and QCD will be part of any realistic calculation and they usually lead 
to infrared singularities. Any method aimed to a numerical evaluation of diagrams must be able 
to handle the infrared problem. For one-loop diagrams we have seen that the method is able to 
extract the infrared pole in dimensional regularization with a residue and a finite part that can 
be treated numerically We have to extend the treatment of infrared divergences to multi-loop 
and we will start the discussion with one specific example: we dehne the on-shell derivative of a 
two-point function, where possibly some of the internal masses are zero, as the d/dp^ derivative 
evaluated at the mass shell of one of the non-zero internal masses. This paper concerns infrared 
conhgurations only and the general case will be presented together with vertices |K , 


Consider now the on-shell derivative of Indeed, for the simplest topology the on-shell 
derivative is infrared finite, see Appendix A. 

It is important to recall that a necessary condition for the presence of infrared divergences is 
that the Landau equations are fulfilled. In the case of we see that s = (mi -|- m 2 ± m^)'^ 
and m 3 = (mi -|- 1712 )^ are satisfied by m 2 = mi = 0, mi = m 3 = m and s = m?. However 
itself is not infrared divergent but self-energy diagrams enter into the calculation of wave function 
renormalization factor and, therefore, also their derivative with respect to is needed. This 
derivative, as in the one-loop case, shows an infrared pole when computed on-shell. Therefore, we 
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consider the case m 2 = m^ = 0 and mi = m 3 = m, a typical example of which is shown in Fig. p. 
For these valnes of the internal masses we get 



Figure 6: A two-loop diagram contribution to the VF-boson self-energy. 


5121 =-2 


X 




F (e) / dx dy dz 
71 I Jo Jo Jo 


121 


x(l — x) (1 — 


'{ml + zY) U;{-> + 
where the quadratic form is 


t /121 = -p z + {p - m )z -b m (1 - )y + 


X 


m 


X 


(119) 


( 120 ) 


First we take the derivative of with respect to p^ and next we go on-shell by putting p^ = —m^ 

5.9.1 On-shell derivative of method I 
Proceeding as described above we obtain the derivative as 


q121 _ 


\ r _„y/2 


m^ V TT m^ 


F (e) / dx dy dz (r(l — x) (1 — y) 


X 


JO 

{z{l - z) [(1 + e) {z^ - x-2-^ + I -f z (3 ^ - 2) |os, 

where we have introduced 


Ao.:12i = + ^ = 


X 


x'-' ■■ l-x’ 

We introduce new variables y = 1 — y' and 2 ; = 1 — 2 ;', so that 


T Tos 2 2 i 2 i y \ 

^121 Ws;121 


and consider the term 


(1 - - ;] (x„a2i) 


-2-e 


( 121 ) 

( 122 ) 

(123) 

(124) 
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To go further, we can use the following relation: 


V, {z^ + = (1 + 6) d, {z^ + f 

V=(dy-^d,,d,, dy), d=[z^-^,-2z, 


After integration by parts we use 


{dz] dy) {z‘^ + ^ = e[-2x-,-^){z^ + ^) 


^ A , y \-l-e 


In this way we derive 


' =-^ 


2 / /i 


\7rm 


[ dx xil — x) ^ S, 
Jo L J 


(125) 


(126) 


(127) 


(128) 


where iS is a combination of one and two dimensional integrals. Due to the presence of an infrared 
pole in the on-shell derivative we must keep e 7^ 0. Next we change variable, y = zy', and introduce 
an auxiliary integral 

Tf .. f J.. f J u f ^ , y \—€ /ion\ 


Using the relation 


we are able to derive 


f dy f dzy^ z'' {z + ' 

Jo Jo X 

1-1 —{ydy + zdz) {z + ^) = 0, 


(129) 


(130) 


/(/i, z/) = 


^2 -\- jjj -\- y — 6 


[ dv[y‘‘0 + ^)-’ + v‘'b + ^)-‘\ ( 131 ) 


In this way the evaluation of S is simply reduced to the calculation of a one-dimensional integral 
which can be cast into the following form: 


S= dySy, 


( 132 ) 


with an integrand of the form 


^y~[ A 


i„-W 2 , 1^1 




+ ^ ( 1 /I -/2 - y2-U2) (1 + X y)-^-^ + i (1 + f )■'• 

First we consider the integrals that produce a pole at e = 0, therefore 


( 133 ) 


[ dyy^'^^'^ ^{l + ay) " = — e [ dy q LA ^ a = X,X (134) 

Jo e Jo y ^ y 
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Successively we multiply the result by [a;(l — and expand around e = 0. The remaining 
integrals can be computed in terms of a long list of master integrals that will be reported in 
Appendix B. The ^-integration gives 


[ dx x{l — x) ^ S = S-i e ^ Ao -|- S'! e, (135) 

Jo L J 

with coefficients 

S _1 = 1 , So = 2 , Si = j-ic(2). (136) 

To obtain our hnal result we have used the following relation: 

Li 2 = LB (x) + lna;ln(l - x) - ^ In^ x - ((2), (137) 

and the following integrals: 

f dx{^^ —Li 2 (a;) ; ln2(l-a;); ln(l-a;)} = {-C(2); C(2)-l; 2; -l}. (138) 

0 cc 

The hnal result for the on-shell derivative reads as follows: 

s™ = ^ _ 2 (2 - t) 1 - 7 + 7 (4 - 7) - 1 C(2) + O (o)]. (139) 

5.9.2 On-shell derivative of method II 


An alternative derivation of the on-shell derivative of will also be presented. We can use the 
following relation 

> + - b 1 (140) 


L l + e 

and derive an alternative result 


= -(z‘^ - 

os ;121 ^ 2;^"^os;121’ 


121 r (1 + < 7 ) / \ 


qLZL _ 

Up 


\7rm 


) lo L lo ~ ~ ^ ~ ^os!m(141) 


Instead of changing variables we can use another important identity, 




(142) 


Once this relation is applied and the integration by parts performed we introdnce new variables, 
y = I — y' and 2 ; = 1 — 2 ;', so that 


Cl 21 _ C" 

^ 2 m^ VTrm^ 


(e) [ dx\x{l-x)\ ^ {SX-S 2 ), 
Jo 


im 
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Si= [ dy{l-y)y ^ {y + X') % 82 = j dy f dzy^/‘^ ^ {z^ + X'y) ", 

Jo Jo Jy 


( 144 ) 


where we used X' = X ^ = 1/x — 1. Next we change variable, y = zy', and utilize the auxiliary 
integral /(/x, z/) of Eq.(|131|), which gives for 82 the following result: 


5o = 


1 — e JO 


dy 9 J (1 + A" 9 ) ‘ + 9 i (y + X')' 


(145) 


From well-known properties of the hypergeometric function we get 


81 =--X^FAe, ] 1 - - ] -X) - 

1 ^ 2 LK , 2’ ^ 2 -e 


^ X"2Fi(6,1-|;2-|;-X) 


x^ 2 Fiie, x) 


*”€(!-€) 2 ’ ^^ 2 ’ “■’^'+(l-e)( 2 -€) 


.YSn(£,l-|;2-|;-A') 


2 e 2 

— -r X^2Fx(€ , 1 : 1 : 1 — x) -j- 7-77-r 

e(l-e) ^ ’ 2 ’ ^ (l-e)( 2 -e) 


x^2Fi(€, x) 


(146) 


After introducing the hypergeometric series, we are able to perform the x integration explicitly. 
Finally we expand around e = 0 by using: 

n—1 ^ 

r (n + m) = r (n) |l -|- - — j) m + O | (147) 

j=i J 

and arrive at the result of Eq. (|139|) . To obtain Eq . (11391) in method II we have used the following 
series: 

00-1 oo-in-i 00-1 

- = c(2), y y ^,, = c(2)-i. (148) 

;S n(n+l) 3 .Sn^(n + 1) 


5.9.3 Another on-shell derivative of 

There is another conhguration of interest shown in Fig. 

b) m 3 = 0, and rrii = m, Vi 3, (149) 

which corresponds to the non-leading Landau singularity cki = 02 = 0 and s = m?. We obtain 


nl 


nl-y 


Sl^ = j r (6 + 1) dx I dy dz 
where we have introduced 


a;(l-a:)l ^ z {1 - z) U \ Ux,y, z), (150) 


os ;121 




m 


x(l — x) 


(151) 


The evaluation of will be postponed till Section where it will be discussed with other 
similar functions. 


27 





























w 


Figure 7; A two-loop diagram contribution to the VF-boson self-energy. 


5.10 Tensor integrals of the family 

Let us define the following function: 


,4 c:1-*/2,2,1-j72 




(/ii... /ij I i/i... lyj) = [ (Tqi d^q2 ^ 

J (gi + mi) ... 


0.1 (I 2 " 


<l 2 


(152) 


{{^2 + vY +'mlY 

where we assume i + j < 3 and where a and fd have been changed in order to account for the 
effective degree of convergence. For tensor integrals we always indicate explicitly the total number 
of internal lines. In any realistic calculation the hrst step is about simplifying numerators and 
denominators as much as possible. For example we will use results as the following one: 

,2e r 


TT^ 


d^qi d^q2 


p-q 2 


{qi + mj) ... ((g 2 + pY + ^l) 

121 




S' (p , mi, m2,7724) -b (m3 - ml -p^) S'4 (p , mi,m2,m3,m4) 


(153) 


After possible simplihcations few irreducible integrals, of the type appearing in Eg. (P-52|) , will 
remain. Their structure will be as follows: 


^l-i/2,2,l-i/2^ 




nl 


i dx [' dy r dz 
0 Jo Jo 


X 


'(p.. .p |p.. .p) = ( — 1 F (e 
\7r s 

n=—l 

The coefficients Q are: 

QYi = -2 xz {yl - z^) QY = xz{l- f) QYi = 0, 
Q^_\ = -2z{yl-zY QY = z{l-l) Q°\ = 0, 


(154) 


and also 


[- 2 a;V-f 

„2 


Qf = (pi - zY 


Q7i = 


e-2 

,^- 2 ^.( 1 -.) 

L e e(e - 2) 

3 5 1 


^Y-^)ipl-zY 


2 e- 1' 


) x^^ — x{l — x) 


+ z^ 

4] 

~ 2 P 


/I 2 

1 - - a;2-4-^-^ , 

e e(e-2) J 


(155) 
























( 156 ) 


= -2 (^2 _ ^ 2 ^^ g 02 ^ ^2 g _ ^ (^^2 _ ^ 2 ^ 1 ^ 

Ql^ = (1 - -) + x (/i^ - 2 :^) -, g^\ = -a; (- - I ^-). 

Note that all terms proportional to X)~ 2 l~^ Q-i contain a factor /r^ — z^ which, in turn, 

allows us to raise the negative power through Eq. (|49|) . 

Several of the properties that we have shown for hold for a more general class of multi-loop 
diagrams, with N >2. 


6 The topology 


is a special example of a general class of two-loop diagrams with iV-|-2 internal lines which 
are overall convergent for N > 2 and contain the aj logarithmically divergent sub-diagram. We 
have 

^4q1ni ^ ^2e J ^ ((gi - ^ 2 )^ + Y[ ((^2 + + m^+ 3 ) \ (157) 

2 = 0 

where the momenta ki are linear combinations of the external momenta pj, ko = 0 and ki = 
Pi + ... + Pi. We introduce Feynman parameters Zi for the propagators of the q 2 loop. 


JV-l 


dz^ = Y[ dzi 6 {1 - Zl) , 


2=0 


iV -1 

Zl=^ Zi 
i=0 


( 158 ) 


and also a {z^j-dependent momentum and mass 


iv—1 


iv—1 


G=Y. =Y, Zi [kj + , 


1=0 


i=0 


( 159 ) 


and obtain 


^ijvi _ _2 ( /^ 

\ vr 


r (iV - 1 + e) 


d[P] [{ml -F x^~ 


+ - 


4-e 


2 Ar -2 + e 


X' 


2—N—e 


,(160) 


where x and y parameters have been introduced to combine all propagators arising after the q 2 
integration. Additional quantities are as follows: 


X = -Py +{M -mjy + m^, 


1 - 6/2 


= / dx I dy I dz^ x{l - x) "'y^ ^{l-yY^. 


(161) 
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Next we change variables, Zi = z'Jy, and obtain 


QlNl _ _2 (/L 

\ vr 


X 


r r dy r- 

e Jo Jo Jo 

l-N-e I ^ ^ ^ ^2—N—e 


N—1 


dzi a;(l 


Xl 


-e/2 


i=0 


{ml - z’Hz) + X 


X 


2 iV-2 + e 

where the qnadratic form x has been rewritten as 


X = z^Hz + 2 K^z + L, Hij = -ki-kj, ^ (^i + "^i+s) > L 


[I - yy/'^ 5{y - Zl) 

(162) 

ml{l-y). (163) 


7 The 5*5 family 

There are two diagrams in the family, given in Fig. ^ and 5^^^, given in Fig. |^. 
is overall logarithmically divergent and contains the logarithmically divergent snb-diagram ay. 

and all its snb-diagrams are convergent. Their evalnation will be discnssed in the following 
sections. 

7.1 topology 

The topology, explicitly shown in Fig. H, is given by 



Figure 8: The two-loop diagram of Eg.(164). 


4 r-131 _ ,, 
TT — fJj 


2e 


d'^qid'^q2 


{ql + ml) ( {qi - gs)^ + mfj {ql + ml) ((g 2 + pf + mfj (g| m, 

If m 3 7 ^ ms then is the difference of two diagrams, 

1 


-. (164) 


5131 = 


mt — mo 


S' (p ;mi,m2,m3,m4) - S' (p ; mi, m 2 , ms, m4) 


(165) 


otherwise it is a special case of with /cq = ^2 = 0, and ki = —p. The qnadratic Xi^i will be 


^131 = 2 :^ - (1 - h4 + + (pI - pl)y + py 


(166) 
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where we have used Z 2 = y — zi — zo and, moreover, z = zi. Limits of integration are 0 < zq < y — zi 
and 0 < Zi < y. Therefore, for m 3 = we obtain 


Cl31 



r(2 + £) 

e 





dz 


a;(l 


X 



^131 


2 -e /2 
1 + e 


X 


-l-e 
131 . 


5 



(1 -9)'''^(!/ - 2 ) 


(167) 


The ultraviolet singularities can be easily derived by taking the derivative with respect to mg of 
5^^^, therefore we obtain no double-pole, as expected, and a single pole given by 


-2 (--A, 


d 


dm% 


21 

/ dx dylnx (x, 1,1/) = — (--A, 
Jo Jo s e 


n l-y 

dy- 


10 


X 


(168) 


121 


Note that the residue of the single pole is what we expect in order to compensate the overlap¬ 
ping divergence coming from the associated subtraction diagram of Fig. ^ The integral can be 
computed, giving 


[' dy- 
lo X 


y 




121 

y± 


h43 


2 

1 

2 L' 


h34 ± A^/2(l,/i^,/i^)l, 


ln(l- 

y+' 


-i„(i--)]}, 


(169) 


where A —A -|- and 6 0+. If we introduce ym = I/ 2 /I 34 and y± = 1 — l/y±, then it follows 

that 


In n - In y- = In F + i,(y_, F), in F = Ay"(l, i4) d _ 

^ ^ ym yVm 


+ 0(A). (170) 


Therefore, for A —a 0+(0_) we observe a square root singularity in the imaginary(real) part. 
Recalling that the counter-term, in the MS-scheme, is — 2 i 7 r^/e and that the remaining integral 
is a Co-function with pi = —p 2 = p we obtain the following result for the subtraction diagram of 

Fig. 


5^31 = 


se '-Jo 


dyXjL + 0(e) 


X 


(171) 


121 


7.2 Landau equations for when mg = m 3 

The Landau equations for this topology and with mg = mg are as follows: 

ai{ql + ml) = 0 , 02 ((gi - ^2)^ + m^) = 0 , 

asiq^+ml) = 0, a4{{q2 + pf + ml) = 0, 

og {ql + mg) = 0, 


(172) 
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Figure 9: The subtraction diagram, containing a one-loop counter-term (represented by a x) associated 
with the two-loop diagram of Fig. 

and also 


aigi/, + a 2 (gi - ^ 2 )^ = 0 , 

—a2{qi — q2)ii + ot^q2ii + O-a (?2 + p)/i + *^5 ?2/i = 0. (173) 

The leading Landan singnlarity occurs for at 7^ 0,Vi We multiply the two equations Eg . (11731) 
by q2fj, and respectively. This gives an homogeneous system of six equations. If all a, are 
different from zero, the singularity will occur for 

9 9 9 9 

= -mi qi = -m|, 

qi- q2 = \ {rrtl -m\- m|) p ■ q2 = \ {s + ml - mf). 

The equations become as follows: 

—2 m\ ai — {ml + ml — ml) a 2 = 0 , 

{ml — ml -|- ml) ai + {ml — ml — ml) a 2 = 0, 

2 p- qiai + { 2 p- qi - s - ml + ml) 02 = 0, 

{ml + ml — ml) 02 — {ml — ml -|- ml) (03 -|- 0:5) -|- (2p ■ gi — ml + ml — ml) 0:4 = 0, 

{ml — ml — ml) a2 — 2 ml {a^ -f 0 : 5 ) + (s — ml — ml) 0:4 = 0 , 

{- 2 p ■ qi + s + ml - ml) a2 + {s + ml - ml) (0:3 -f- 0:5) - {s - ml + ml) 0:4 = 0. (174) 


Compatibility between the hrst three equations requires the conditions 


mt 


ml = (mi m2) , P- qi = — 5 - 5 - 7 {s + ml- mj). 

mf — ml + m 3 


(175) 


As a consequence, it follows that 


m 2 i 

ai = -q;2i 0^2 = - 7 -;-r 

mi m 2 (mi -|- m 2 ) 


1 r 


|-(mi m2)^ (0:3 + 0:5) -h I [s - (mi m2)^ - m 


0(4 


}■ (176) 


is a solution, if and only if s = (mi -|- m 2 ± m 4 )^ and m 3 = mi -f- m 2 , therefore, representing the 
leading Landau singularity. We are now in a position to attempt the evaluation of this topology. 
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7.3 Evaluation of for = m 3 

In Eq.(|167|) we have a power —2 — e but for the hrst iteration we prefer to use the identity 


1 + ^131 ' = ^ 


-2-e 
131 ’ 


and obtain 


= — ( — ) r (1 + e) / dx I dy I dz x{l — x) ' {1 — y) 


1 / p 


nl /-l 


l+e/2 


e/2-1 


S \ TT S 


/O JO 


X (y - \x,y,z), 

with X = x{l — x) X- It is more convenient to change variables, y = 1 — y' and write 


^131^ 


1 / /i 


S \ TTS 


2 \ 


r -1 ^1 


r (1 + e) / dx dy 


nl-J/ 


dz 


10 Jo 


a;(l — x) 


l+e/2 


y 


e/ 2-1 


X {l-y-z)X^^^ ^{x,l-y,z), 


where now the quadratic form is 


and where Eq.(|^) can be used. We introduce 

^131 "^(l! h'3’A^4)' 

The hnite part of becomes 

Q131 =-^ y- J131 

fin oh ■> 

* '^131 n=0 


/ 3 I 31 = r dx [' dy [' ' dz 
Jo Jo Jo 


3 \nX {x,y,z) + (2 z - 1 - - ^ 3 ^) 


1 . 2 N 




y 


m = 


rl ^1 


dx / dy {in (x, 0 , 1 /) [2 (1 - ?/) + (2 1 / - 1 - - yl^) (x, y) 


10 Jo 


2 

^■34 


/InX (x,i/, 0 )' 


1 - n 


y 


lnXi 3 i(x,i/, 0 )]}, 


= - P 34 j dx ln(x — x^) — InX^gj^ (x, 0 , 0 ) InX^g^^ (x, 0 , 0 ), 

fl31 _ ^ 

k - 5 . 

where L.,, is given by 


LiX^,y) = lii(l -y)+ lnx + lii(l -x)- In Vai {x, 0 , y) 


131 


(177) 


(178) 


(179) 


''■'131 = “*)lii 3 i> X,3,(x,l-y,z)^ z - (1 + ft-1*4)2 + (1*1 - 1*3)!/+ 1*3. (180) 


(181) 

(182) 


(183) 


(184) 
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The obvious advantage of this result is that does not depend on x. However, = 0 for 
s = (m 3 ± 1 x 14 )“^. In this case we have a quadratic form + bz + cy + d where Ad = b‘^. Therefore, 
we consider hrst the case ^34 = —2 /i 3 ; we have 


Xi3i ( 2 ^. V, «) = (« + + (A - A) V- 


(186) 


It is convenient to change variable, z = z' — so that the integral to be evaluated becomes 
^ lo lo (1 - 2/ - ^) y, z) 


= / dy 




'M3 


dzy‘1^ ' (1 + jiz - V - z) (z^ + cy) ‘ 


(186) 


with c = — yl- After changing variables, z^ = t we use 


(t + cy) ^ = - -dt{t + cy) ^ 


(187) 


and integrate by parts. The hnal result, which is valid for ^/s = — m 3 and 7714 > m 3 can be 

cast into the following form: 


J =-- 


In /i 


2(1 + yz) 


’•l+fJ.S 


e +3 “'M3 

1 /■! . /■i+/'3-M dz r 


dy 


Iny 




2 Jo 


dy 


'M3 


ln (^2 + cy) - (1 + ys) 


ln( 2 ;^ + cy) 


y 


(1 + + 3 ) 

1 + 


'■ 1 +M 3 rln?/ln(l +/i3 — ?/) hi^?/-| 

dy 


'M3 


2/i3 Jo 


dy + cy) - 


y y 

|^ ln(/ii + c|/) ^ 


+ 


1 2 2 

In 


2/^3 


(188) 


The case ^/s = m 3 — m 4 , with m 3 > m 4 is obtained after the substitution z = 1 — z' in Eq.( p. 86 | ). 
We obtain 


<7 = ^ dy dzy^/"^ ^ {z -y)x^l^\x,y,l- z) 

/*! /*1+A^4 

= dy dzy^/'^-^ {z - 114 -y){z^ + cy)~^~% 

t/ 0 J y-\-fj,4 


'y+M4 

+131 (^ 12 / 1 1 - ^) = ^^ - ^ + {yl -yl)y + y’i 


(189) 


and the case considered here corresponds to P 43 = —2 1 x 4 . Therefore the pseudo-threshold s = 
(m 3 — m 4 )^ is covered. 

Always for = m 3 we have a second method which is based on the observation that the 
following identity holds: 


+131 


/ dx{l — x) 
Jo 


d 


dMl 


1) 5'33(a:V; mi, m 2 , My,), 


(190) 


34 




















with /i^ = Ml/s and x + The idea is to change variable x —> /i^ in the integral 

and perform integration by parts. Unless nnmerical differentiation is used, this method will fail 
at = 0 and, therefore, it will be no longer discussed. 


7.3.1 Evaluation around s = (m 3 + 

In this region we use 

/^3 ~ ^ (h'34 ~ = -^(1) h-S’^'4); (191) 

and write the master integral J of Eq. (|186| ) as 

J= dy dz {1-Zm-y - z){z‘^ + cy - - (192) 

JO J-Zm 4 

where Zm = and c = As a first step we will set A = 0 whenever possible; J is the 

sum of three contributions, the first one being 

pi rl — Zm—y 

Ji = — dy dz z {z"^ + cy)~^~^. (193) 

0 'J—Zm 


Here we use the identity 


- d^{z^ + cy) ^ = 2z {z^ + cy) 

e 


-l-e 


which gives, after integration by parts, the following result: 


A = -ln-^ 
e 1 - z„ 


+ 2 hi^ Zm - ln^(l - Zm) 


f\yl^'^v{y,-Zm) lny{y,l-Zm-y)] 

2 Jo 


y 


y 


where ri{y, z) = z^ + cy. The second contribution, J 2 is dehned by 


J, = 


1 fl-Zm-y 

dy / dzy^/^z^ + cy)-^-\ 


JO J Zm 

Here we will us another identity, i.e. 

--^dy {z^ + cy)~^ = c{z^ + cy)~^~\ 
giving, after integration by parts, the following result: 


(194) 


(195) 


(196) 


(197) 


J 2 


1 

c 



dz 


lnr]{l - Zn 


z, z) — In 77 ( 0 , z) 


(198) 
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The last contribution is the most difficult to obtain. We start from the dehnition of J 3 , 

rl pl — Zm—y 

J3 = - dy + cy)~^~\ (199) 

JO J—Zm 

The y integration range is further divided into two regions, [0 ,1 — Zm\ and [1 —Zm, !]• In the former 
one 1 — Zm — 2 / is positive and we write 


rl — Zm rli 

J3i= dyy^/‘^-^'^ dz{z^ + cyy 

Jo Jo 

where W = z^ and Y 2 = I — Zm — y- Therefore J 31 is a combination of new integrals denoted by 

Hi= f dz {z"^ + cy)-^-^ = {cy)-^-^Yi 2 ^ 1(1 + e, ^ ^ ( 201 ) 

JO 2 2 cy 


2 fYi 


\-l-e 


( 200 ) 


Using a well known property of the hypergeometric functions |16| we obtain 




y 


13 Y-^ ( r 

^ ’ 2 ’ 2 ’ cy’ 

f c \ 13 c 

f 172 y ) 2-Ui(l + e, - + e; 71 + ^) ~ 172 ^)’ 


yi-j " ’ 2 ’ 2 ’ Yy 

where the coefficients B 12 are as follows: 

r(3/2)r(-l/2-e) _^ 

^ r(l/2)r(l/2-e) l + 2e’ 

where 'ip{z) is the Euler psi-function [^. In this way the functions Hi become 

Hi = BiYi ^ 2 Fi{ 1 + e, 2 +*^’ 2 ~*”*'’ ~ ^ 


( 202 ) 


(203) 


(204) 


At this point we observe that the term proportional to B 2 will give rise to a divergency after the 
?/-integration. If A is kept different from zero then Eq.(p04|) becomes 


^?i = Sil'r'-"‘2n(l + e, 1 + e; 5 + e; -+ ^2(cj/- 1(205) 

Consider now the ^/-integration for the second term: we obtain 


f' 1 Z-n 


H,2 = 


dyy^d^ ’{cy--\) " 
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= 2 




A 


2 -Fi(- + e, - ; 1 + 2 ’ ~ 

A 


^ 21 (1 - .„)-V^-./^ .-./.- 1 , , 1 ^ ^ ^ __, 


(206) 


where the coefficients fei ,2 are 


, r(-i/2-£/2)£ , r(i + £/2)r(i/2 + £/2) 

"1 — Tn 1 ^ -TTvT 7T) "2 — - 


r (1/2 -e/2) 2’ r(l/2 + e) 

For the term proportional to hi there is no UV pole and we may set e = 0. Therefore we get 


(207) 


H,2 = -2{l - ^ 


) + 2 ^c-/^ 


= -4A-i/2 


arcsm 


'2 ’ 2 ’ 2 ’ 4c(l-z^) 

ucd^l2 <- [7 -1> - ^ - ^( 7 )] ■ <™) 


Setting A = 0 whenever possible produces the following result: 

Hi 2 = - 2(1 - ZraY^^"^ + 2 (- A - i 5)“^/^ ^ - In c - ln(- ^ * 5) - 7 - '0(^) , (209) 


where 6 0+. Returning to Eq.( 205 ) we consider the part proportional to Bi and the corre¬ 

sponding ^/-integral which becomes 


/* 1 — ^7n / 1 

= fi{e,y) fi{e,y) = 2 Fiil + e ,- +e ;- +e y). 

We split the integral into two parts, according to the following recipe: 

Ha = 


( 210 ) 


<iy/.(..0) + y-'d—Id 


= -(1 - Zraf/"^ fi{e, 0 ) + 


y 


y 


+ o{e) 


( 211 ) 


The function /j is dehned by 


/,(e, 0 ) = Fr'“'^( 2 / = 0 ), /,( 0 , 0 ) = 


Y^{y = 0y 


ffn \ 1 P 71 1 3 cy 1 1 / 2 , Yi + {-cyyF 

/.(0,,) = - 2 Fi(l,-;-;-^) = -(-q/) 


( 212 ) 
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Since J 31 is the following combination, 




(213) 


2=1 


we obtain the hnal result by further inserting 


Bi ..- 1 - 


2^Yr^-^^{y = 0){l-Zmr^ = 


- y = - — 


2 rl 


'-e 


H“ ^ ln(^l ^m) 2 In ZjYi 2 




9 r1 S 

^- - - ln(l - - 2 


The second part of the function is J 32 dehned by 

J32 = f [ dz{z^ + cy) 

Jl — Zm J— Y 2 


-1-e 


(214) 


(215) 


where we can set e = 0 from the beginning. We obtain 


J?.2 — - 


1 /■! dy 


1 3 




c J\ — Zm y j — 




2=1 


cy 




(216) 


In all formulas the square roots must be understood with the replacement cy ^ cy — id, where 
5 ^ 0 +. 

Let us return to Eq. (|204| ) and consider also the term proportional to Bi for A 7 ^ 0. In this case 
we will use again Eq. (pllj) but with 


1 


cy - A/4, 


— " 2 Fi{l + e , - +e; - +e; - 

Due to the pole at e = 0 we must expand /j(e, 0) to O (e). We derive 


(217) 




13 r(3/2 + €) 

2 n(i+f. 2 +<; 2 +<;!/.)-r(i /2 + j) 

where we have introduced yf = A/( 4 Y'j^). Therefore, expanding around e = 0 , we obtain 

+ e F, 


^ / 1 3 n, 1 , 1 + Xj 

2n(i + f. 2 + <; 2 + ';!'<) = ^‘‘'rvv 


F=- 


dx - 


2 Inx — ln(l — 


yfx^ - 1 


, 1 , I+ Xi 

H-In --, 

Xi I- Xi 


(219) 
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a result which is based on the standard integral representation of the hypergeometric function [|l 6 |. 
The integral can be computed exactly; introducing Xi = dyf = \fX/{2Yi) we obtain 


1 r. 


- 

X, 

1 


2xi 

1 

2xi 


Lia {xi) - Lis (-Xi) 
In 2 ln(l + Xi) + Lis 
In 2 ln(l — Xi) + Lis 


4xi 
I + Xi 
2 

1-Xi 


ln^(l + Xi) — ln^(l — Xi) 


Lio hr 


- Lis hr 


, 1 , l + Xi 

H-In --, 

Xi 1 - Xi 


( 220 ) 


where Lis(1/2) = (C(2) — lh 2 )/ 2 . To obtain Eq. (p 20 D we have used the following integrals [|: 


dx ■ 


Inx 


_i — Lis (T Xi) , 


rl 


X ± X. 


ln(l ± Xix) 


dx 


ln(l ± Xix) 


= ± 


/ dx ■ 
lo It XiX 


1 r 

= T — 

Xi 


In 2 ln(l T t) — Li 


1 ± XiX 2 X 

It X. 


hh(l ± Xi), 


+ Lis ( - 


( 221 ) 


( 222 ) 


Furthermore, Lis ( 1 / 2 ) = (C(2) — lh 2 )/ 2 . The implementation of this result requires distorting 
the a:-integration contour for Js (Eq.(|198|)) in order to avoid the points 


x± = 


2 /i| L' 


L? - L2 + Ls ± L?) L 2 ) 


(223) 


These two points coincide for /i| = (/xi T /i 2 )^- In practical cases this may happens for mi = m 3 
and ms = 0 where we have 


,.2 _ Li „ _ ,.2 ,.2 _ ,.2 I ^ 

) c Lx Ls Li 

X X 


(224) 


Actually there is no problem at a; = l(c = 0) since the result for Js is 


J 2 


1 X 
fj^i 1 — X 


/ dz |ln h — In 

— Zm. 


2 I 2 

^ + Li 


1 — a; 

X 


(1 - Zm 



(225) 


Our result for the threshold region deviates from the main approach of expressing the integrand as 
products of polynomials and logarithms of polynomials and requires some analytical work before 
attempting the numerical evaluation. However, in this way, we have been able to isolate the 
square root singularity which shows up in the coefficient of the single pole and will disappear 
after including the corresponding subtraction diagram. The net result can be safely computed by 
numerical integration. 

^ We disagree with Eq.(l) p. 34 of ref. ^ 
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7.4 Derivative of and infrared poles 


Self-energy diagrams enter into the calculation of wave function renormalization factor and, 
therefore, also their derivative with respect to is needed. We concentrate on the on-shell 
derivative of where = m 3 and compute 


d 


-S'^ (p ;mi,m2,m3,m4,m3) = 


131 _ cl31^„2. 

dp^ 

{p^ + m\- ml) mi, m2, m3, m4, m3, m4) 


2p 

+ S'^^^(p^;mi,m2,m3,m4,m3) - mi, m2, m3, m4, m4) 


The result for the derivative is 


2 \ £ 


551 = - 


TT s / e Jo 


dx 


a;(l — x) 


-e/2 


dy 

10 JO 


X 


{r {3 + e)z{y- z) {pi -pl-l) [{pi - ^ 2 d 


2 + e 


2-^/2 i d 


+ r (2 + c) (p - 2 .) ^ }. 


1 + e '"131 


(226) 


(227) 


It is important to recall once more that a necessary condition for the presence of infrared diver¬ 
gences is that the Landau equations are fulhlled. In the case of we see that s = (mi-|-m 2 ±m 4 )^ 
and ml = (mi -|- m 2 )^ are satished by the QED-like conhguration, m 2 = m^ = 0, mi = m 3 = m 
and s = m^ but, for instance, they are not when mi = m 3 = m, m 2 = 0 and s = m^ but m^ = M. 
This conhguration, that corresponds to electron (or W-boson) self-energy with 2 = 7 and 4 = Z, 
is therefore free from infrared poles. 

Let us consider a realistic case. Fig. which represents one of the contributions to the electron 
self-energy in QED. In this case the conhguration is 



Figure 10: An F^^^-type electron self-energy in QED. 

(228) 

(229) 
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1 ) mi = m 3 = m and m 2 = m^ = 0 


and, therefore 


^131 = ^ 


{1 +p^)z+{p^ - pl)y +pI, with 


2 _ 

X 















(230) 


with = m?/s. The derivative is needed for on-shell electrons, therefore /i^ = 1 and 


Tos;131 ~ ^ ^ x' 


We therefore obtain 


1 1 

to 


\ 7rm^ / 


r(2 + e) /•! 


xfe-2^) [(;-x=)xW,+ 


dx 


2:(1 — x) 


-e/2 


dy [ dz{l- 


^os;131 




(231) 


1 -|- £ ''os;131. 

The well-known infrared singularity shows up in parametric space at y = z = 1. Note that the 
diagram itself is not infrared divergent but only the on-shell derivative, due to the factor y — z in 
the integral of Eg. ( |167|) . To continue our derivation, we use the identity 


1 + 


1 


1 + e 

to write the on-shell derivative as 


{ydy -h zdY 


X 


-l-e 

os;131 


= 1 --. )X 


- 2 -. 

^os;13l’ 


(232) 


5. 


1311 

P \ OS 




rrr \ vr 


r( 6 ) 


dx 


x{l — x) 


1-^/2 


dy [ dzil-yY^^ 


X (y — 2 z) 3-I---I- {ydy -|- zdz) 




-l-€ 

os;131 


(233) 


After integration by parts we obtain 


5, 


1311 

P \ OS 




2 V TT 


r(i + £) 


dx 


x{l — x) 


1 - 6/2 


y. j\lv f\lz(v-2z)(l-vY'^-\-J-;^ 




2 \ TT m? 


r(i + £) 


dx 


a;(l — x) 


- 6/2 


X 


dy [ dz{2z-y-l) ^ + (- l)y 


1 - 1-6 


The total result is the sum of three contributions. We write 


J= dy dz{2z-y-l)y^/‘^ ^ z"^+ {--l)y 

Jo Jy X 

and introduce X = 1/x — 1. For Ji we obtain 

Ji = - - f dy j dz dz {z^ + X y) 

€ Jo Jy 


X 


1 - 1-6 


— Jl — J2 — Jo 


(234) 


(235) 


1 P 


e Jo 


dy \yJ^-^ (1 + Xy)-^ - y-J^-^ {y + X)" 


(236) 
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For computing the remaining two integrals we use the following relation 


1 + (- zd^ + ydy) (z^ + Xy) 


\-l-E 


1 + e '2 

After integration by parts we obtain the following results: 

nl 
10 


= 0 . 


(237) 


{l-e)J2= f dyy/^l + Xy)-^-^ + y-^/^y + X)-^- 


rl 


10 


(1 + 6) Js = - / dy //2-1 (1 + X y)-^-^ + 


The Ji integrals are successively folded to become 

ri 


Ii= dx 
Jo 


x(l — x) 


1-^/2 


Ji. 


We illustrate in details one particular case and refer to Appendix C for further details: 


J= / dyy 
Jo 


-e/2-l 


(.y + X) 


-e-1 


This integral gives rise to an hypergeometric function [K 


J=-- X-^-^ 2FAl + e, 

e i u -r , 2 ’ 2 X^ 


(238) 


(239) 


(240) 


(241) 


Using well-known properties of the hypergeometric functions we obtain 


J = -- x^^^ 2 Fi {1 -f e, 1; 1 - ; a:) 

2-^1 (1 + ^5 


l+e 


=- X 

e 

+ Aa (1-x)-^-^/2%Fi(-^e, -^e; l-x) 

= [Ai 2 Fi(l + e, 1; 2 + |e; 1 - x) + A 2 x^^^ {1 - x)-^-F^^ 


where the coefficients A, are given by 

r(l-e/2)r(-l-3/2e) 


Ai — 


We next compute 


r(-3/2£)r(-£/2) 


2 + 3e 


Ao. — 


r(l-e/2)r(l + 3/2e) 

r(i + e) 


X = / dx 
Jo 


a;(l — x) 


- 6/2 


iX — Xi -|- Xa- 


(242) 


(243) 


(244) 
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The first integral is 


2r(l-6/2)r(l + 3/2e) /•! 

' 6 r(i + 6) 


[ dxx^^’^^l — x) 
Jo 


-l-2e 


2 r(l-6/2)r(l + 3/2e)r(2 + 6)r(-26) 2^^ 

e r(l + e)r(2-e) e ^ ' 


(245) 


In expanding around e = 0 we use 

r(-2e)^-^-7 


72 + C(2)]e + 0(e2), 


r (a + A) = r (a) {l + ^/>(a) A + i [ij\a) + V’'(a)] + O (a") }. 


It follows that 


12 11 ^ 
^1 — “o “I-1" 2 + -r C(2)- 


The second integral will be 


2 .^2, T (2 + 3/26) T (n + 1 + e) 

2 + 36 ^_g T (1 + 6) T (?7, + 2 + 3/26) 


+ +7++^ ^h': + ++' £ (1 - 


(246) 

(247) 

(248) 


where we have introduced the hypergeometric series. After integrating over x we can put 6 = 0 
and obtain that 

J2 = 1-C(2). (249) 

We also have the subtraction diagram of Fig. ^ The counter-term, in the MS-scheme, is 
After taking the derivative with respect to we put = —m? and obtain 




131 


2 

( ] 


\ vrm^ / 


r(2 + 6/2) + 


/ dxx{l — x) 
Jo 


-2-e 


(250) 


In both cases one has to approach the mass-shell in a interval of the dimension n where the integral 
is convergent and continue the result in the number of space-time dimensions. For instance, from 
Eq. (p50|) we obtain 


il31 


1 

to 






r2 , ,1 1 1 9 1 ^ ^ 

^-(1+7)- + 1 + 57+j7" +jC(2) + 0(e) 


(251) 


Consider now another conhguration, depicted in Fig. 

b) rrii = m 2 = m, m 3 = = 0 and = M. (252) 

It is immediately seen that b) does not correspond to the leading Landau singularity but rather 
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Figure 11: A second S'^^^-type electron self-energy in QED. The innermost loop can be a fermion of any 
flavor. 


to the sub-leading one corresponding to ai = 0:2 = 0 and s = We start from 

^ .2\ n i^y 

10 JO JO 


= Jo dy ^(y-"{x,y,z), 

(253) 


take the derivative and put obtaining 

= fwl r(2 + e) [ dx [ dy dz\x{l-x)] (y-z) z{l-z)U-^-1{x,y, z), 


TT 


10 Jo Jo 


os;131 


where the quadratic form U is 


u^s-ioA^^y^^) = ^ ^+^x{^-y)^ ^x = 


m 


os;131 

Now we make use of the identity 

zil — z) (x, y, z) = 

and integrate by parts obtaining 

1 


x(l — x) 


- --- - (1 — z)dzU ''(x,y,z) 


(254) 

(255) 

(256) 


Cl31 _ 

^ nnh 


r (1 -|- e) / dx x{l — x) ^ {Sh -|- 5c) 
Jo 




(m 


2 \-l-e 


2M^ \7i 

where we have introduced 

Sb = {ml)-^-^ f dyy{l-y) 

Jo e[2 + e) 

‘5c = ^ dy dz{l-yy/^-^{2z-l-y)Uy^yy{x,y,z). 

In order to compute Sc, we make the transformation y = 1 — y' and introduce the function 

r-l rl-y 


(257) 


(258) 


E{y,iy,a) = J dyj dzy^^ z^ Uy^^^^^{x,y, z) 


(259) 
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which can be computed explicitly with the help of the identity 




(260) 


giving the following result: 


I' dyy^ (l-yni + y) 17“ ^ 3 ^ (x, y,l-y) (261) 


2 , 0 . 2 ^ ly 3 Jo 

Substituting this result into Eq. (|257|) and expanding around e = 0 we obtain 


9131 ( h 

2M4 [ttM^ 


r (1 + e) / dx dy S _2 e ^ e ^ + 5( 


(262) 


where the coefficients S are given by 

5_2 = -4, 5_i = 2X 

[3 + 2/ o2/ 


■y_ _ 

f+i 1 

^x 

\yX^y 


+ 


3/^2 


5n = 7f 


X 


2 + L + 2(—) +2L+1 -4 + C(2) 

V \7/ A// + J 


X 


^ 2 /X 




3/^2 '3 


(263) 


where we have introduced x .131 = X (1 — y)"^ + ^2 _ ^ 


L = lnx--lnX-- Inj/, L+=fi^ 

2 2 + Ux^ 


2 V 2 /X, 


(InX + In y). 


(264) 


Note that we can treat these coefficients numerically since x is a semi-positive dehned quadratic 
form with zeros at y = 1 and x = 0 or x = 1 which are compensated by the numerators (note 
that we have changed again variable, from y to 1 — y). A convenient trick consists in extracting 
these factors in order to obtain a manifestly hnite expression. Consider, for instance, the following 
integral, 

/•I ^ 2 ^^ _ ^^2 


X = / dx dy ■ 
Jo Jo 

after some straightforward manipulation we obtain 


X" 


(265) 


X = / dx dy 
Jo Jo 


3 rx'^y{l — xyY x(l — x) x'^y{l — xy)- 


Xi 


+ 


xi 


+ 


xi 


X1 =x(l-x|/)(l-?/)^ + p^ 


, n 2 


X2 = 

X3 = a: 


l-xy{l-x) + y y, 


l-y{l-xy) +/i^ 


(266) 


where it is immediately seen that the three quadratic forms are strictly positive over the [ 0 , 1 ]^ 
region. 
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By looking at Eg . (11501) we see that the evaluation of is very similar to the one we have 
just performed. Using the same technique we obtain 




nl 


-e/2 


= — 2 — r (e) / dx dy dz a;(l — x) (1 — y) 


.e/2 


TT 


10 JO 


X 




where the quadratic form now is 


U = -p z + (p + m )z + (1 - y) 


(267) 


(268) 


First we take the derivative with respect to p^ and next we impose the on-shell condition, p^ = 
—obtaining 


0121 _ 
^pb 




nl /•! 


^l-U 2 


- 


\7im^ 


r (e) / dx dy dz a;(l — x) (1 — y) 
Jo Jo Jo ^ -I 


X {z(l - [(1 + e) - .Y) x-/-' + 2 x-y] + -- (3 X - 2) x-‘y}l< 

where the on-shell quadratic form is 


Xos;121 =^' + ^(1-2/)’ ^ = 


a;(l — x) 


To proceed, we can use the following relation: 


,-2-e 


1 + ^ (yd, + zdO\ y „ x 2 i ■ 


(269) 


(270) 


(271) 


In this way we derive 
r (e -|- 1) / 


q121 _ 
^pb 


r-1 ^1 


dx I dy r dz x{l-x) (1 - p)"/^ ^ z(l - z) y j " (272) 




j Jo Jo Jo 

Now we can use 

■r M 

^os;12l7 2e ^ ^'-os;121 

Once the integration by parts is performed we introduce y' = 1 — y and it follows that 


os;121 


^(1 -^) (Tosn 2 i) ' = - 7 T -(1 


(273) 


(- 4 ) / 

^ 2,m^ Xirm^ I Jo 


1 nl 


dx I dy\x{l-x)\ ^ iyXy"-Xo"+ I ^ dzx " )(274) 


1 - 2 / 


where the various y’s are expressed as 

X = z‘^ + Xy, Xy = {.^-yf + ^y, Xo = Xy. 


(275) 


46 



























In order to perform the ^-integration we use the function -E(/i, z/, a) introduced for in Eq. (p59|) 
and derive: 

f dy [ ^ dz {1 - x~^ = [ dy {1 + y) Xy" (276) 

Jo Jo 1 — e Jo 

After the integration of the xo term, we introduce Uy = x {1 — x) Xy = x {1 — x) {1 — yY + y and 
perform the e-expansion, obtaining: 


Cl21 

‘^pb 


1 

2m? 



r(i + e) {s^ + sY), 


(277) 


where the coefficients are 


2 r2(l + e/2) 
e2 r (2 + e) 


^2 



dy 2 InUy + 




(278) 


7.5 topology 



Figure 12: The two-loop diagram of Eq 


The evaluation of is more involved than ones already presented for the other two-loop 
self-energies. We have 


TT 


^221 ^ ^ 2 e j ^ ^ 2 ^ ^ ^^2 ^ ^ 2 ^ _ ^^^2 ^ ^ 2 ^ 1 

X (ql + ml)-\{q, + pY + ml)-\ (279) 

First of all the a — j propagators are combined with parameters Xi,X 2 - After performing the 
gi-integration we introduce two additional parameters y, z, perform the g 2 -ffitegration and derive 


^221 _ _ / r(l + e) 


rxi 


dxi I dx2 

Jo 


X2(l - X 2 ) 


1 -l-e /2 


X 


[dy 


^221 


(1 - y)pa, + zp + z{p^ + ml) + {y- z)ml + (1 - y)ml^, 


(280) 
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where xi, 2 -dependent mass and momentum have been introduced: 


Xi - X2 

Px = - p, 

I- X2 

o 1 


"" X2il-X2)^ 

The quadratic form in y, z has coefficients 


{xi - X2){1 - Xi+ X2)p^ (1 - Xi)m\ + {xi - X2)ml + X2ml . (281) 


H = 


-pI p-Px 


P-Px 


-p 


.2 > 


^ (2Px + ml- -2p-p^ + ml- mfj , 


L = m: 


XX t'x 


pI-15. 


Since p^ and p are simply related by 

p^ = Xp, 


X = 


Xi - X2 

I- X2' 


(282) 

(283) 

(284) 


the matrix H is singular. To overcome this problem we introduce new variables y = y', and 
= z' + Xiij' — 1)., giving for the quadratic form the simple expression f /221 = az'^ + bz + cy + d 
with coefficients 


a = —p 


c = ml - ml^ + X(p^ + mg - m|) d = 
After the transformation we will use the identity 


b = p^ + ml — ml, 


m: 


X{p^ + mg — ml) 


l + -d,+ 
e 


" ' 2e 




u-^ = 

221 


B 


221 


4 a 221 


with i? 22 ]^ = 6^ — Aad. This factor reads as follows: 


i?22i = K~P : ^i) + 4:P {mi^ - ml) - 4 Xp {p + m^ - ml) = sb. 


221 ' 


(285) 


(286) 


Before continuing in the evaluation of the diagram we study its singularities. 


7.6 Landau equations for 

The Landau equations corresponding to are 

oi {ql + m\) =0 02 ((<?! -f p)^ -f ml) = 0, 

«3 ((<?!- 92 )^-fJ^i) = 0 ai{ql + ml) = Q, 

«5 {{q2 + p)'^ + mf) = 0 
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and also 


+ Oi2{qi + p)^ + tt 3 {qi — q2)^ — 0 , 

- 03(^1 - ^ 2 )/. + aAq2^l + as (g 2 + = 0 . 


(287) 


The leading Landan singnlarity occnrs for a* 7 ^ 0,V7 We mnltiply the two eqnations Eq. (p87|) 
by ^ 2^1 and respectively. This gives an homogeneons system of six eqnations. If all a, are 
different from zero, the singnlarity will occnr for 

9 9 9 9 

qf = -mi qi = -ml, 

qi- q2 = \ (m.| -m\- m\) p ■ q2 = \ {s + m\ - ml), 

p- qi = \ {s + m\ — ml) 

The eqnations become as follows: 

—2 «! + (s — — m2) a2 — {m-^ + m^ — m^) = 0, 

/2 2 | 2 \ 1 / 2|2 2 \ /2 2 2 \ n 

— — rrtg + 7714) «! + (s — + m3 — mg) 0:2 — (m^ — m3 — m4) 0:3 = U, 

/i 2 2 \ / 2 | 2 \ i /2 2 2 | 2 \ n 

(s + m^ — m2) «! — (s — m4 + m2) a2 + (m4 — m2 — m4 + mg) 0:3 = (J, 

/ 2|2 2 \ /2 2 | 2 \ 1 / 2|2 2 \ n 

(m^ + m3 — m4) 0:3 — (m4 — m3 + m^} 0:4 + (s — m2 + m3 — m^) ctg = U, 

(m^ — m3 — m4) a3 — 2 m4 0:4 + (s — m4 — mg) ag = U, 

— (m^ — m^ — m4 + mg) 0:3 + (s + m^ — mg) 0:4 — (s — m^ + mg) 0:5 = 0. ( 288 ) 


In order to have non trivial solntions the determinants of the coefficients mnst be eqnal to zero. 
In both cases the solntion is given by 


s = 


± 


2 ml 


{2/2 2 2 2 2\ / 2 2\ t 2 2\ 

■jmg (m3 — m^ — m2 — m^ — mg) — (m^ — m4) (mg — m2) 


\ / 2 2 2\ \( 2 2 2\1 

A(m 3 ,m 4 ,m 4 ) A(m 3 ,m 2 ,m 5 )J 


(289) 


which represents the so-called anomalons threshold for this topology, the leading Landan singn¬ 
larity. Actnally the singnlarity occnrs when the prodnct of the two lambda-fnnctions is positive. 


i.e for 

m3 < mQ4 or mQ2 < m\ < m^g or mg > mQ4 ( 290 ) 

where anxiliary masses are defines as 

mQ4 = min{(mi — m4)^, (m2 — mg)^}, ( 291 ) 

mo2 = min{max{(mi — m4)^, (m2 — ms)^}, min{(mi -|- m4)^, (m2 -l- ms)^}}, ( 292 ) 

mQ3 = max{max{(mi — m4)^, (m2 — mg)^}, min{(mi -|- m4)^, (m2 -|- mg)^}}, ( 293 ) 

mQ4 = max{(mi m4)^, (m2 ms)^}, ( 294 ) 


The snb-leading Landan singnlarities follow from a stndy of the rednced diagrams of Fig. 


49 
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4 



+ perm.; 



5 


+ perm. 


Figure 13 : The reduced diagrams corresponding to the topology of Fig. 1 ^. Their leading Landau 
singularities are the sub-leading ones for 


7.7 Evaluation of method I 

To compute we must, first of all, considers the roots of 6221 
that 6221 =0 for 


dehned in Eq.( P86|) . It follows 


xi + (hs + 7 - (1 + h 5 - l^l)XlX 2 - (1 + /ii - nl)xi -t- (1 


■ h2 + hs ■ 


■/i4-- A)x2 + /i? = 0, (295) 


where A = A(l, /ig). Before discussing the general case let us consider some simple example, e.g. 
TTii = m,Wi 3 and three cases: m3 = 0 , m3 = m and also m3 = M, where M is some large mass, 
M > m. The curves corresponding to i? = 0 are given in Fig. 0 where the integration domain, 
0 < Xi < 1 and 0 < X2 < Xi is also shown. Before raising powers we introduce = SX221 > with 


A221 ^ 


h45 Z+ifil- /xL - ^ h«) y + l^lx + ^ h45) 


where, as usual, = 1 + jof — The raising procedure requires 


X, 


-l-e 

221 


^221 


1 + — dy + — {z + Z)dz 
e 2e 


X 


221' 


(296) 


(297) 


where the object of interest is 

^221 ~ ’ ^- 4 ) hs) ~ 4 (hxx ~ 1^4) ~ 4 X p, 45 , Z = — — X = — , ( 298 ) 

/ i — X2 

and 0 < 1/ < 1 , with X{l—y) < z < y+X{l—y). After a Laurent expansion in e we will have terms 
proportional to powers of lnx22ij following these operation we substitute back z = z' — X{y — 1 ), 
so that the order of integration will be specihed as 


1 




dxi- 


(299) 
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Figure 14 ; The xi — X2 plane with the curves 6221 = where 6221 dehned in Eq. (| 286 |) . Here rrii 
m, Vi 7^ 3 and m3 = 0 ,m,M with M > m. 
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The argument of the logarithm is 


X, 


T221 


221 


^221 — B^Xi + C*^, 


^2(1 - X 2 y ’ 

where the coefficients are given by 

= (1 - y) (1 - X 2 ), Bl = -(1 - I/) (1 - X 2 ) (/i54a;2 + fiu), 

C° = X2(l - X2f {z^ - X 2 ) [{y^l - i4)xl - ^^ 2^2 - yl 


(300) 


(301) 


and = 1 + — /i|. Zeros of 6221 avoided by distorting the xi integration contour in a way 

that preserves the correct imaginary part of In ^22^^. As done for the sunset diagram we consider 
the two roots of the equation X22j^ = 0 , say If the roots are complex conjugated we select 

a contour which starts at Xi = X2 and bypasses the zeros of in the upper (or lower) half-plane 
to return to Xi = 1 . 

If the roots are real we must distinguish between two cases: when is positive then the cut is 
between xi^ and xia, while for negative the cut is [—U [xi^, -|-cxo]. Furthermore, with 
Xi = a + i f 3 , the imaginary part of the logarithm is f 3{2 aA^ + B^). This consideration immediately 
tells us the sign of the imaginary part when xi approaches the real axis on the cut. This sign is 
crucial when we distort the contour in the complex plane since, for xi real and on the cut, the i 6 
prescription, gives — in for the imaginary part. 

Let us consider one of the normal thresholds, s = (mi -|- m3 -|- ms)^, as we have done for other 
diagrams. This corresponds to the three-particle cuts of Fig. |^. Let Xi{l,r} be the two points 
where X„^-. = 0 and Xi± the two roots of = 0 . Whenever xi± are on the real axis we want to 


distort the integration contour. 


In Fig. IT we give the real and imaginary parts of these four points as a function of X2 for hxed 
y and different values of 2;. For decreasing X2 the branch points and are complex conjugated 
and no distortion is needed. For values of 2; near to y the branch points become real when xi± are 
still complex and, again, distortion is not needed till the value of X2 where Xi± are real. However, 
for some value of 2;, that we call Zth_{y) the following happens: the two branch points pinch the real 
xi-axis exactly at Xi^ = Xi_ and distortion is not allowed anymore. A part for having a larger 
number of Feynman parameters, the situation is completely similar to the one described for the 
sunset diagram [Q. Since there is no ultraviolet divergency we may set e = 0, obtaining 


1 /*! p 

5*^^^ = — / dxi / dx 2 (1 - X2) dy 

s Jo Jo Jo Jji 


1 rv+^y-y) 


X(l-y) 


dz X ^. 
221 


Our hrst step will be to use the following relation: 


y ^ = 

"^221 


{l-[ydy + ^{z + Z)d,\\Yix22i}- 


(302) 


(303) 


"221 


We obtain 


^221 _ 


4 A /■^i A ry+^y-y) r r 1 1 '1 

— / dxi / dx2{l-X2) / dy dzll- ydy + -{z + Z)dz lny22ir- ( 304 ) 

^221 dxy-y) L 2 j 
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Figure 15 ; The three-particle cut of 5 ^^^. 


where = X 2 (l — 2 : 2 )A. Three new quantities are needed: 

^^21 = X22i(^ = + -2/)) = (1-A:)V- [l-h5 + AL-2X(l-X)]r/ + /i2^ + X^ 


A221 = X22i(^ =-^(1-1/)) = X y + [/i4-/i^^-2X \y + y^^ + X , 
aLi = A221 (2/ = 1> ^ = 1/) = 1/^ - h 45 1/ + 


^221 

Therefore, using 


we obtain 


^221 _ 


rl py+Xil-y) 

dy dz = 

3 JX(l-y) 


pX. p\ r\ p\ 

/ dz I dy+ dz dy, 
Jo 21-4 Jx 


(305) 


(306) 


rl pxi 3 rl ry+Xil-y) n 

/ dxi / dx2\- dy dz In X221 “ 

Jo Jo '■2 Jo Jx(l-v) Jo 


^^221 


Xil-y) 


dy \nx 


221 


- i 1^' dy [[X-y + Xy + Z) InxJ^^ - {X + Xy + Z) Inx^^i] + 


(307) 


Next we change variable, 2 ; —*• z', so that 0 < 2 ;' < 1 and write all x’s as polynomials in xi with 
coefficients that are functions of X 2 ,y and z. In particular we introduce 


= 


^221 


221 2:2(1-2:2)2’ 


Xl,,=A\x\ + B\x, + Cl. 


(308) 


we hnd 


A^ = Y I-X 2 + YX 2 , Bj = -Y {1-X2)yl2 + ‘^Xx2 , 

=Y{1- X 2 ) {-2 xjy + X 2 yl 2 + xl + yj) + Y‘^xl + X 2 {l - 2 : 2 )^(/ig -Y)y 


A001 = 


X2 

221 


221 X 2 (l- 2 : 2 ) 2 ’ 


Xl^^=Alxi + B-^,x, + Cl 


t2_2 


= F (1 - 2:2 + 1 ^ 2 : 2 ), B'^ =-Y (1 - 2 : 2 )(/ii 2 + 22 : 2 ) + 2 ^ 2:3 

c| = F (1 - 2 : 2 )(p. 322:2 + 2:2 + flj) + Y‘^x\ + 2 : 2(1 - X 2 fll\y. 


(309) 

(310) 

(311) 
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Finally we introduce 


X 


0 

221 


_ 221 

0:2(1 -X2Y' 




221 


Apl + Bill + Cl, 


(312) 


where 


= (1 - y) (1 - 2 : 2 ), = -(1 -y){l- X2) {,yliX2 + / iia ), 

Cl = 0:2(1 - 0:2)^ {z^ - + yly) - (1 - i/)(l - 0:2) \{yl - yl)xl - y%X2 - yl 


( 313 ) 


Here we used 

Y = l-y. ( 314 ) 

At this point we must realize that the situation for is different from all previous cases. Indeed 
and all its sub-diagrams are convergent and, for this reason, only logarithms appear after the 
hrst iteration and not logarithms squared. Therefore, the strategy for distorting the xi integration 
contour will be the following. Let us introduce variables Aj = 0:1 — Xu, where Xu is any root of 
the argument of a logarithm, and split the logarithms according to 


In [A^xI + B^xi + 


= InA^-F?7(A^, A^Ajj)-h?7(A^, A«)^ \n{xi-xu), 

i=LyR 


(315) 


where rj is the usual Veltman’s function, 17(0, b) = ln(a6) — In a — In 6, and logarithms are dehned 
with a cut along the negative real axis. The final expression for the diagram will be the sum 
of terms proportional to r^-functions and terms of the form ln(a:i — Xii)/Kii. The former do not 
create a problem in distorting, when needed, the integration contour. For the latter, however, if 
one of the roots of A^, or both, is real and internal to the interval [0,1] then the distortion for the 
corresponding term must be examined with care. If the imaginary part of Xu is positive (negative) 
then we move the contour into the negative (positive) imaginary half-plane, so that the cut of the 
logarithm will never be crossed. If xu is real and negative no problem will arise and we can distort 
in any of the two ways. If its is real and positive then the distortion is hxed according to the i 6 
prescription: for instance the roots of axf + bxi + c — i 6 corresponding to 6^ — 4 ac > 0 are such 
that 

Imxiqi^fl} = =Fsign(a) 5 . ( 316 ) 

The reason why we cannot apply this argument to the other cases discussed before is linked to the 
presence of In^-terms. After the splitting of the logarithm we will encounter terms proportional 
to ln(a:i — Xi^) ln(a:i — Xia). When the roots are complex conjugated we will have cuts both in 
the positive and in the negative imaginary half-planes and the possibility of distorting the contour 
ceases when they approach the real axis. 
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7.8 Evaluation of method II 

An alternative evaluation of follows when we change variable in Eq. (|28U| ), y = 1 — y'^ so 
that 



j•y+X(l—y) 

dy / dz{l — yY''^ {az^ + hz + cy + d)~^~'^ = 

Jx{l-y) 

dy J dzy'^'"^ {az^ + bz — cy + d')~^~'^. 


(317) 


This method has the advantage that d' = c + d = mY with a corresponding change of h —> h' , 
where — d/i^, is xi, 2 -independent. In this way the xi —X 2 integration is free of numerical 

instabilities. The only point in the external parameter space where a failure will occur is for 
s = (m 4 ± m^Y, representing normal and pseudo thresholds associated to the two-particle cut of 
Fig. 0 They are non-leading Landau singularities. The fact that 6^21 i® ^ independent means 
that, with this method, we can perform the x-integration without having to distort the integration 
contour. As observed the method fails at the non-leading Landau singularities corresponding to 
the one-loop complementary sub-diagram of (in this case there are two sub-diagrams with 
the same number of lines). When we are in this kinematical conhguration and s Y (^1 ±^2)^ the 
strategy will be to interchange the order of the operations, q\ ^ q 2 - When (mi±m 2 )^ = (m 4 ±m 5 )^ 
the method fails. Using this approach, we set e = 0 and get 


e22i 


X 


1 rxi 

- / dxi / dx 
s Jo Jo 

/-I 
"221 


dy 


rl+{X-l)y 

Ixy 


dz 


X2(l-X2) xJAxi,X2,l-y,z- Xy), 


(318) 


X221 X2,l-y,z- Xy) = - yl^z+ (/iL - yl + X ylY y + yl. 

Successively we use the identity 


X, 


-l-e 

221 




1 + - a + 


2 e 


a 


X 


221 


A(l,/il,/ii) u ' g 

The derivation continues, as usual, with an integration by parts, where we use 

rx pzix .(z-i)/{x-i) 

/ dz / 

10 JO Jx Jo 

Next we dehne 


rl |^l+{X-l)y rX rz/X i-l r{ 

dy dz= dz dy+ dz 

Jo JXy Jo Jo Jx Jo 


dy. 


2 _ 
t^xx 


I'Z 


0:2(1 - 0:2)’ 


(319) 


(320) 


(321) 


(322) 
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and introduce four quadratic forms: 


X221 = ^2(3^1 - X 2){1 - X 2 )ynl 5 + 2^2(1 - X 2){1 - Xi)y{yl^ - 2 ) 

+ 0:2(1 - a:2)^(l - yyl + yl- yl^) + 3:2(1 - xify"^ + (1 - X2)yul^, 

X221 = ^2(3:1 - X2fy‘^ + 3:2(1 - 3 : 2 )V 4(1 -y) + {l- X2)yul^, 

X221 = ^2(3:1 - 3 : 2)(1 - X2)yy% + 3:2(1 - 3 : 2)(1 - Xi)z{yl^ - 2 ) 

+ 0:2(1 - 3 : 2)^(1 - yyl + yl- yl^) + 0:2(1 - o:i)^2:^ + (1 - X2)yi^l^, 

X221 = ^2(3:1 - 3 : 2)(1 - X2)y%iy - z)+ X2iXi - X2fz‘^ 

+ 3:2(1 - 3 : 2 )V 4 (i -1/) + (i - X2)yjyl^. 


Note that, in the limit 0:2 1, we obtain 


1 2/1 \2 2 3 4/1 \2 2 

= (1 - 3 : 1 ) y , x,,i = = (1 - 3 : 1 ) ^ 


"221 "^221 


221 "^221 


Similarly, in the limit X 2 0, we obtain y. The result for is as follows: 


^221 _ 


where the function is 


s\{l,yl,yl) 


nl rxi I 

/ dxi / do :2 — J 221 (^i’^ 2 ,l/, 2 :), 

00 00 0:2 


J 221 Jo L '^221 Jo '^221 ’ 


and 


1^1 yjl-xi) p 

221 2 (1 - 0 : 2)2 


1 ^^ 221 -1^4i 


^ ^45) In X221 + l(y + l ^^45) In x' 


221 


2 _ 3 1 - 0:1 

"'221 


lnxj 2 i -Inx: 


221 


2 (1 - 0 : 2)2 L 

It is straightforward to show that in the limit 0:2 = 0 


1-0:2 L2 


I Inx^ 


221 


+ 1 


so that we obtain 

o221 


J22^{xi,0,y,z) = 0 , 


4 

sX{l,yl,yl) 



1^ /221 (xi, 0 : 2 ) ^ 


(323) 

(324) 

(325) 

(326) 

(327) 

(328) 

(329) 
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The integrand is defined as follows: 

/?■ = 

with auxiliary functions given by 




(330) 


P = 
221 


X2 


3 

2 

3 

2 L 
1 


1 


(^1’ ^2, y, z) - In (a:i, 0, y, z) 
[ (1 -x^f ^221 (^1’ ^2, y, z) - In {xi, 0, y, z) 
lnX22i(^i’^2,l/,^) -lnX22i(^i>0>l/>^) 


4i = ^1/(1 -^i) [(TZ^ \Tix\^^{x,,X2,y)-\nx\ 


1 


j^yp-xi) 

+ X (y + mL) 


x52idi>dy) 

1 

^221 (^1 > ^2, y) - In d y) 

^221 ^2, y) - In X221 dn 0, y)' 


( 331 ) 



Figure 17 : One of the two-particle cuts of 5 ^ 2 ^. 


7 . 8.1 The case s = {m^ ± ms)^ 

As already pointed out the derivation fails in this case because P45 = ±2 ^4. The alternative 
solution is quite similar to the one adopted for other topologies. If P45 = — 2 p 4 we can change 
variable z' = z + so that 

X 221 ^ + (yL - ^4 - 2 X/i4) y, (332) 

with Z- < z < z+ and 


z- = X y + ij. 4 ^, 


z+ — 1 + (X — T)y + fX4, 


( 333 ) 




























where z± >0. We change again variable = t and nse 

{t + ay)~^ = dt\n{t + ay), (334) 

and integrate by parts obtaining 

1 f^+ 1 /"i r 1 

- / d?/ / dt\n{t + ay) = - / dy — \n{z\ + ay) 

I Jo Jzl I Jo '-z^ 

1 1 

-ln( 2 ;^ + ay) + - j dt ln(t + ay) . (335) 

If /i 45 = 2 p 4 this derivation is of little nse becanse the new limits of integration can be negative. 
However, in both cases we can write 


P-4 — ^ h45 '^(1)^45 hs) 


(336) 


introdnce Zm = ^ 45/2 and transform z' = z — Zm- With z. = X y — Zm, and 2 :+ = l + {X — l)y — z„ 
we obtain 


which gives 



(337) 

f dyR-^ [in^+^-ln^^], 

Jo ^ Z- — R Z- + R^ 

(338) 


where we have introdnced 


= a = hL-h4-2W/i2, 


(339) 


showing that no singnlarity will appear for A = A(l, /i|, pg) = 0. 


7.9 Derivative of and infrared poles 

In this section we stndy another infrared divergent object, the on-shell derivative of the 
integral corresponding to the conhgnration a) where 

a) m 2 = m 4 = 0, mi = m 3 = mg = m and = —m3. (340) 

It is important to recall that a necessary condition for the presence of infrared divergences is that 
the Landan eqnations are fnlhlled. In the case of we see that the system of Landan equations, 
Eq.( |288D is satished for this conhgnration of parameters if 03 = —Oi and 05 = Oi. itself is 
not infrared divergent but the on-shell derivative is. Let us consider the diagram of Fig. p! 8 | that 
corresponds to conhgnration a). Solid lines correspond to the particle of mass m, while dashed 
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lines give the massless particle. After taking the derivative with respect to we will have two 
contributions where the propagators (gi — + rn^ and (g 2 + p)^ + will appear with power —2 

respectively. To study qualitatively the infrared behavior we neglect qf, and qi ■ q 2 compared 
with p ■ qi and p ■ ^2 in propagator denominators, the so-called eikonal approximation. We further 
dehne a superhcial degree of infrared divergence for an arbitrary scalar two-loop diagram 

= 2n-l-2nm - Um, (341) 

where is the number of massless lines and the number of massive lines. For the diagram 
of Fig. |l^ (after p^-differentiation and for on-shell p) we obtain pijn{Sp^^) = —1. This shows the 
presence of an infrared pole singularity. We have 



Figure 18: A special case of topology. Solid lines denote a massive particle, dashed lines a massless 
one. 


,221 _ C ' ^221 _ I \ f f 


Sp = ITT = - - r (2 + e) / dxi 


dp 


TT 


2:2(1 -2:2) 


1 -l-e /2 


rl rl + {X-l)y 

x dy dzy^!^ 

Jo JXy 


z{z -1)- py 


221 ’ 


where the quadratic form is dehned by 

^221 = ^ 2 /) + ^Ix 

Xi - X 2 


X = 


1 - 2:1 

P = -A, 


1 - X 2 ’ X 2 

and is given in Eq. (|281|) . Next we dehne the on-shell limit of 1 / 221 ’ 

TT ! ‘2 2 \ 2 2 I (1 — 2:1 -f 2:2)^ 

^221 (p =-m)=m Xo,. 221 ’ Xos ;221 +vy, y= • 

Inserting this result into Eq. ( |342|) we obtain the on-shell limit of the derivative, expressed as 

,2 1 /■! 


(342) 

(343) 

(344) 


(345) 


S, 


221 ;os _ 




rxi 


7rm^ 


T{2 + e)— dxi dx 2 0 : 2(1 - 2 : 2 ) S 


Jo 


-l-e/2 


pi pl+{X-l)y 

^ 221 = / dy / dzy^' z{z-l)-py 

JO JXy *- 


{z^ + py) 


-2-t 


(346) 
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In the evaluation we must keep e 7 ^ 0 and one cannot expand around e = 0 after the y — z 
integration because of the presence of overlapping infrared divergences; here we face a novel 
feature and the procedure is made even more complicated. Our solution is derived by adapting 
the general algorithm of [^. Before explaining in detail how the residues of the infrared poles 
and the hnite part can be cast into a form which allows for a reliable numerical integration we 
notice that the 2 ;-integration can be performed in Eg. (|346|) after using the identity: 


1 


^iy9y + ^d^)] + vvT = 0 - 


(347) 


We dehne the following integrals: 


hj = dy 


rl+{X-l)y 

Ixy 


dz z^y^ {z^ + yy) 


-2-e 


(348) 


and obtain 


I 20 — 


Iio — — 


1 — e Jo 
1 


loi — 


e Jo 
1 


dy Q-^-^ {i + iy- iV - eV) + (1 + O' 

dy Q-^-^ (1 - iV) + y-^'" L-^-^ (1 + O' 


dy Q-^-^ (1 - iy) + y-^'^ (1 + 0 


L = y + X^y. 


1 — e Jo 

with ^ = X — 1 and where 

Q = i^ + ^yf + vy, 

The relevant combination of integrals is 

dJ = I20 ~ dio ~ pdoi- 
Collecting the various terms we obtain 

I ' dy {-^ [2/^/2 (1 - e'2/') + 2 /"^/' L -^-^ 


(349) 

(350) 

(351) 


c _ 
*-^221 


+ 


1 -e 


X 


y"" {i + ((-p)v-{(-p) (y" - (V) + y-'" {X y-p)x l-^-] }. (352) 


From its dehnition, Eg. ( |345|) , we see that y = 0 for xi = 1 and X 2 = 0. While the guadratic form in 
y contains a constant term, the linear form y + X'^y vanishes if the parameters {x'l = 1 — xi,X 2 i 2 /} 
are set to zero. This is the origin of the overlapping divergence that we must disentangle. First 
of all we change variable in order to factorize the divergent terms: Xi = 1 — x'^ and,successively, 
X 2 = Xix' 2 -, so that divergences will arise from the behavior of the integrands at the point x) = 0 . 
After the transformation we get 




3^1 

1 - 0 : 10 : 2 ’ 


y = xi 


{I + X2? 

0:2(1 - 0:1X2) ’ 


(353) 
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When we exchange the xi — X 2 integrations the limits become 0 < 0:2 < + cxo and 0 < xi < 
1/(1 + X 2 ). We change again variable, Xi = x '^^/(1 + X 2 ), so that 0 < < 1 and 


e = - 


1 + a;2 - X1X2 ' 


7] = Xi 


{l+X2y 

X2{1 + X2- X 1 X 2 ) 


(354) 


7 . 9.1 Defoliation of IR divergent integrals 

At this point we introduce a simple example to explain the method. Consider 

/ = / dx [ Q = x{l + x) + y'^. ( 355 ) 

Jo Jo 

A minimal set of zeros is a minimal set of parameter values that satisfy Q = 0. In our simple 
case we have {x = y = f)}. A defoliation is the procedure that allows for a factorization of the 
divergent part (in the limit e —0) through different steps. In our case, 

I=[l dx f dy+ f dy [ dx} = f dx f dy [x~^~^ + y~^~^ ( 356 ) 

where the defoliated Q’s are: 

Q{x=o} = ^ + x + xy‘^, Q{y=o}=x{l + xy)+y. ( 357 ) 

Since the hrst form never vanishes over the whole parametric space, including the border, the 
procedure stops here for the hrst integral, while for the second one has to iterate until the condition 
is met. To give an example, if we start with Q = x{l —x) +y‘^ the hrst defoliation gives x{l —xy) +y 
but also l — x+xy"^ which has a set of zeros at {x = 1, ?/ = 0}; then one has to transform, x = 1 —x', 
and iterate. 

Note that the algorithm is perfectly dehned for all integrands that may vanish only at the 
boundaries of the parametric space. is the hrst example of a strategy which dihers considerably 
from the minimal BT approach. In the latter the main attempt is toward raising powers in the 
integrand while here the powers are untouched and the whole idea is about factorizing the singular 
behavior into simple factors that can be integrated analytically, leaving non singular terms to be 
integrated numerically. 

Other cases, where the Feynman integrand is not semi-positive dehned (presence of thresholds) 
are more difficult to analyze. However this never happens for two-point functions and their 
derivatives. Infrared singularities are hnally extracted by using Taylor expansion, typically 

[ dx [ dyx~^~^ f{x,y,e) =-- [ dyf{0,y,e)+ [ dx [ dy ■ ( 358 ) 

JO Jo e Jo Jo Jo ^ ^ + 

Note that only now one can expand around e = 0. Various branches in the procedure are not 
uniquely dehned and one is led to the search of the path which brings the minimal number of 
iterations. 
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7 . 9.2 Sp‘^^ method I 

After these preliminar considerations we go back to Eq.( |352|) . After all transformations we rewrite 
Q and L obtaining 

Q = {1 + X 2 - xiX 2 y a;2 + (1 + xl) {I + X 2 - X 1 X 2 ) xiy + x\ X 2 y"^, (359) 

L = (1 + X 2 Y (1 + a;2 - X 1 X 2 ) Xl + (1 + X2)^ (1 - xi)^ X2 y. (360) 


The integration range is from 0 to 1 for xi and y and from 0 to cx) for X2. It is easily seen that Q 
and L may vanish only for zero valnes of the parameters. In order to extract the divergences we 
nse the relation: 


rx rl ry 

dx dyf{x,y)= dx dy f{x,y) + dy dxf{x,y) (361) 

Jo Jo Jo Jo Jo Jo 

= dx dy {x f{x,xy) + y f{xy,y)) (362) 

Jo Jo 


Since onr fnnctions f = Q,L are polynomials to a power —2 — e, we have /(x, xy) = x~‘^~^g{x, y), 
where g has a zero at x = 0 of lower order. We split the integration over X2 in two parts, the 
first in the interval [0,1] and the second in [1, cxo]. In the former we iterate the procednre nntil we 
reach polynomials that do not vanish inside the parametric space. In the latter we perform the 
procednre only for xi and y. 

In the limit X2 —> cxd some of the integrals may lead to a divergency becanse of the presence of 
divergent factors, other than in Q and L. In that case we perform the transformation X2 —»■ I/X2 
and iterate the procednre for X2 G [0,1]. At the end we obtain integrands where the divergences 
are contained only in the following factors: 


Xi 


,-l-e 


Xi 


„-l-2e 


Xn 


,-1-e 


or y 


-l-2e 


(363) 


The hnal answer is obtained by using the following relation: 


/ dyy^^^f{y)=f dyy y f{0)- 

Jo Jo ^ y ' + 2e 

In this way, after the expansion around e = 0 we obtain 




dX2 


X2(l - X2) 


-1 -l-e/2 


5, 


221 



e 


221 


+ 5, 


221 


(364) 


(365) 


where , 5^ and 5° are expressed in terms of smooth integrands. Furthermore, is the 
221’ 221 221 ^ ^ ’ 221 

sum of integrals which can be calculated analytically: 


52 

221 


1 

'2 Jo 


dxo 


10 


dy 


X2 


1 

2 


roo pi 

dx2J dxi 


1 + X2 
X2 


(1 + X2 + X2y) 


-2 


n 


I + X2 


(xi + X2 + X1X2) 


-2 
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1 /-I /•! 1 I 

-- dxi dx2— -(1 + a;2 + a;ia;2)"^ - - / dxi (1 + xi)“^ 

/JO JO 1 + X2 /JO 

1 /■! 

(1 + a ;2 + xia; 2 )~^ - 77 / rfx 2 (1 + 0 : 2 ) 

/ JO 


^ / dxi I dx2-^ 

2 Jo Jo 1 + X 2 

^ f dxi ( dx2 - 

2 Jo Jo 1 + 0 : 1 X 2 

^(In2-l)-^(ln2-ln^) + ^(ln3-21n2-^) 

- (31n2 - 21n3 - -) = -1. 

2 2 


-2 


( 366 ) 


For the single IR pole we obtain a relatively long expression that will be reported in Appendix D. 
The expression of 5°^^ is too long and of little interest from the analytical point of view. Therefore 
it will not be reported here. 


7.9.3 


method II 


To check the correctness of method I we have also used a different strategy. In Eq. (|346|) after 
changing variables as described after Eq. (P52| ) we exchange the order of integration, so that X 2 G 
[0, cx)] and Xi G [0,1/(1 + X2)]. Two other changes of variables will bring all integrals in the interval 
[0,1]: they are Xi = x)/(l + X2) and X 2 = l/x^ — 1. After rationalization we end up with two 
polynomials to the power —2 — e. They are 


Q = yxi{l - X1X2) [x2 + (1 - X2) J + y xiX2{l - X2) + (1 - X1X2) x 2 , 

Ld = yxlx2{l - X2Y + yx2{l - X1X2Y + xi(l - X1X2) [2x2(1 - X2)(l -y) +xl + {l- X2) 


(367) 


We have written the two polynomials in a form that makes clear their semi-positive dehniteness 
over the cube 0 < {xi,X 2,|/} < 1. For more than two variables, xi.. .x„, defoliation is based on 
inserting the following identity, 

n n 

1 = H n ^ (368) 

1=1 j^l=l 

in the integral. Several defoliations are needed and we summarize the results in the following list. 
First Q*, 


{Xi = X2 = 0} 
{Xi = X2 = 1} 

{y = X2 = 0} 
{y = xi = 0} 
{y = X2 = 0} 
[y = Xi=X2 = 1} 
{xi = X2 = 1} 


Q* —> Qi + 

Ql Qll + Ql2y 
Qii ^ Qiii A Q1121 

Qi2 Q12I A Qi22j 
Q2 Q21 A Q22) 

Q21 ^ Q211 A Q212 A Q2131 
Q22 Q22I A Q222-I 
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( 369 ) 


and then Lh 


{y = X2 = 0} 
{y = xi = 0} 

[y = Xi = 0} 

{xi = X 2 = 0} 

{y = Xi = X2 = 1} 
{y = Xi = X2 = 1} 
{xi = X2 = l} 


Q221 

Q222 


L\ 


-^11 

ft 

^12 

- 

^2 


Q22II Q22\2i 
Q222I Q2222J 

' t I ft 


Ll + Ll, 


L\, +L* 


12 ) 

t 


11 

-^111 + -^112 + -^ 113 ) 

01 + -hi 90 + L\ 


4i + L 


22 - 


(370) 

In the above procednre it is nnderstood that the appropriate change of variables is performed to 
bring the zeros at the origin at each step. At the end of any chain all polynomials are strictly 
positive dehned in the nnit cnbe. One hnal step is needed to disentangle integrals of the following 
form: 

J{a,b)=[ dx [ dy{l-x)~^~^^{l-xy)~^~’'^f{x,y), (371) 

Jo Jo 

where f{x,y) is analytic in the integration domain. We obtain 

J{a,b) = [ dx [ dy {1 + y — xy)~^~’^'' f{l — x,l — xy) 

Jo Jo L 

+ (1 + a; - xy)-^-’^^ /(I -xy,l-y). 

Finally all divergent integrals factorize, an example being as follows: 

K{a,b) = [ dx [ dy x~^~°‘''y~^~^''F{x,y), 

Jo Jo 


(372) 

(373) 


A-(aa) = ^ f (0.0) - 1 / d, (^)^ - i / * ( 

F{x,y) 


F{x, 0) 


X 


l+ae / _|_ 


dx t 
lo Jo 


dy(- 


q^\-\- 0,€ fy\-\-be / _|-1_ ^ 

where we have introdnced a '+' and a ' + +' distribntion 




Jo X 


1 — ae\nx + O 


(374) 


(375) 


'■ t iy ( ) = r f dy Mx.y)-M 0 ,i/)-Mx, 0 ) + /i( 0 . 0 ) 

1 Jo Va;l+ae|/l+6 e/++ Jq Jq xy 

X 1 — e (alnx + 61 ni/) + (9 


(376) 


Also in this case the residne of the donble IR pole can be compnted analytically and gives —1. 
The explicit expressions for the residne of the single IR pole, as well as for the hnite part will not 
be reported. 
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7.9.4 Another configuration for 
There is another configuration where Eq. 


is satisfied, 


b) mi = 1712 = = ms = m, m3 = 0 and s = 0. 


( 377 ) 


Actually this configuration gives a nice example of a necessary condition which, however, is not 
sufficient. Let us consider the corresponding derivative. 


d 


C^^i|6 _ " c221|6 _ 

- 


TT 


2\ e 


r (2 + e) / dxi 


rxi 


X 


dy 


rl+{X-l)y 

Ixy 


dzy’"^'^ z{z — 1) 


py 


u 


dxo 


-2-e 
221\h ’ 


X2(l -X 2 ) 


-l-e/2 


where the quadratic form is defined by 


^22116 ^ -p z + p z + (m^^ - m - Xp^)y + m^. 


and m^^ is given in Eq. (|281|) . Next we define the —> 0 limit of U, 


22116 ’ 


(378) 


(379) 


m 




^221j6|os’ ^221|6|os ^2) y + ^ 2 . 


(380) 


Although the point p = 0:2 = 0 is the candidate for generating an infrared divergence we can 
prove, by direct calculation, that the result is infrared finite. When we set e = 0 it is easy to see 
that the p, z integrations are trivial and 


‘522116 — 


rx\ 


7 ! 7 nr T ^n(^l) Bn{xi) ^ 


(381) 


where the various coefficients are: 

Ae = 3 — 6x1 + 3x1, As = —— + 14xi —A4 = 7 — —Xi + 2^1’ 

A3 = -2 + ^xi, A2 = ^, Bt = 3- 6x1 + 3x4, 

Bq = —10 + 17xi — 7 x 4 , i?s = 11 — 13xi + 3 X 4 i ?4 = —5 + 3xi, B^ = 1 . 


(382) 


The non trivial part of the X2-integration is performed recursively by using the following relation: 

Xi 


J —n. 


/ 

Jo 


, In X2 

dx2 


(1-X2)’^’ 

^ J—{n+l) (^2. 1) J—n T 


J_2 = ln(l - Xi) 
1 1 


n — 1 n — 1(1— xi) 


“ + 


- Inxi. 

Xi — 1 

Xi Inxi 

il-xiY' 


(383) 
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Finally the xi-integration uses the following results 

rl /•! 1 1 

/ dxx~‘^ ln(l — x) = — / dx{x~^ + -x~‘^ + -x~^) 

Jo Jo 2 S 

rl /■! 1 3 

/ dxx~^ ln(l — x) = — / dx{x~‘^ + 

Jo JO 2 4 

/ dxx~‘^ ln(l — x) = — / dxx~^ — 1, 

Jo Jo 

[ dxx~^ ln(l — x) = —C(2), 

1 


11 

18’ 


rl I rl 

/ dxx^lnx = - -/ dxx"^ \ii(l — x)) = 

Jo in +1)^ JO 


/o ' {n + lY' 

The hnal result is 

'’22116- ^2’ 
proving that the b-conhguration is infrared hnite. 


1 1 

n + l^i I 


<SooA \h — ~ 


13 


(384) 

(385) 

(386) 


7.10 Tensor integrals of the family 

The scalar integral has been computed at n = 4 and the absence of double-logarithms has 
shown some advantage for the integration in the complex plane. On the contrary, tensor integrals 
of the same family are ultraviolet divergent and this feature seems lost. However, this problem is 
encompassed by the procedure of scalarization. Scalarization can never be complete, in contrast 
with the one-loop case [^, and we will stop the procedure whenever an irreducible numerator 
is encountered. A tensor integral will be irreducible when the numerator cannot be expressed 
as a linear combination of those propagators which are present at that particular moment of the 
reduction. In a sense we scalarize the obvious and explicitly compute the rest. 

Let us dehne the following function: 


TT 


4 ^2 */2,i,2 j/2^^^ J d'^q2 


Ml M? Ml 

Ql ■■■ Ql <l2 ■■■ ?2 


(gf + mi)... ((g2 + vY + mf) ’ 


(387) 


where we assume i + j < 3 and where a and (3 have been changed in order to account for the 
effective degree of convergence: for i + j > 2 the function is overall divergent, for i > 2 the ay 
sub-diagram is divergent, etc. Note that we always indicate the total number of internal lines. 
Let us consider some example, 

_ (388) 

^For the first equation we disagree with the result of p. 20, ref. |3^ 
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and perform the reduction. Since there are enough propagators, the reduction is particularly 
simple, 


qi-p= 2 {[(^1 + "^2 

— 

' 2 1 2 
Qi + 

q2-p= ^ {[{q2 + p)^ + ml 

— 

' 2 1 2 

q 2 +mi 


^ rr^ ^2 I 21 r 2 I 21 r 2 , 21 2 , 2 , 21 

<?1 ■ <?2 = - ^ I [(gi - q2) + - [qi + - [^2 + -rn^-\-m^ + 7714 j. 


( 389 ) 


etc. and give the following result: 


8 5 , 


1/2,1,2 


(p, p, p I 0 ) = {ml — m2— p y S (0 I 0 ) 


+ {mj - ml- p^Y I 0 |p, mi, m3, mi, m^ 

+ 4 p^ (m^ — m^) — Tp"^ — (m^ — m^)^ | 0 | — p, m2, m3, ms, mi) 

+ 2 {m\ - ml- pY sY‘^’‘^’^{p | 0 | p, mi, m3, mi, ms) 

+ 2 ( 5 p^ - mi + m^) | 0 | - p, m2, m3, ms, m 4 ) 

5°^^(p, p I 0 I p, mi, m3, m4, ms) - (p, P | 0 | - p, m2, m3, ms, m4) 


( 390 ) 


Another interesting example of scalarization is represented by 

8 55 ’^’^/^(p,p I p) = -P 12 P 45 I 0) + p^2P45 ^(0 I 0 I P) m3, m4, mg) 

- P 45 iPi2 + 2p^) 5]^^(0 I 0 I -P,m2,m3,ms,m4) 


P12 

2^45 


^4 50 I 0 Ip, mi,m3,mi,m2) - 55^50 I 0 I - p,m3,m3,m2,mi) 


5y^’^’^(p I 0 I - p, m2, m3, ms, mi) - | 0 | p, mi, m 3 , mi, mg) 


5/2,2.1/ 


+ (p + m2 + P12) 53 (0 I 0 I -p,m2,m3,m4)+P12 [53 (0 | 0 | p, mi, m3, mg) 

— 53^^(0 I 0 I 0, mi, m3, m4) + (p^ — m^) 53^^(0 | 0 | —p, m2, m3, mi) 

- (pi2 + 2p^) 53^^(0 I 0 I 0,m2,m3,ms) -2 53'^^’^’^(p | 0 | p, mi, m3, mg) 

+ *S'3'^^’^’^(p I 0 I - p, m2, m3, m4) 

( 391 ) 

where we have introduced the following combinations: p^j = p^—m?+m^. Other relevant examples 
can be obtained by permutation of masses and all of them indicate that we need at most the scalar 
5s-function plus a combination of tensor integrals with a number of internal legs which is < 4. 


7.11 A realistic example 

To discuss some realistic example of tensorial reduction in the 5s-family let us consider one of 
the two-loop diagrams that in QED contribute to the photon self-energy, see Fig. We derive 


























Figure 19: A contribution to the photon self-energy in QED 


the following expression for the vacuum polarization: 




( 392 ) 


For the function S we simplify the numerator and obtain 


S'= 4 (e — 2)p^ + 8m^ + m, m, 0, m, m) 

+ 2 8m^ — (8 — 2e — BQ{p^]rn,rn) + 8 {2 — e) Bl{p‘^]m,rn) 

— 16 2m^ + (e —2)p^ m, 0, m, m) + 4 e (2 — e) m, 0, m) 

+ 8 (e — 2) S'^^^(0; m, 0, m) + 8 (2 — e) Ao(m) B^ip^] m, m). 


(393) 


Note that S is ultraviolet divergent and must be associated with the corresponding subtraction 
diagrams of Fig. 0. 



Figure 20: The one-loop subtraction diagrams, containing a one-loop counter-term (represented by a x), 
associated to the two-loop diagram of Fig. 19. 


8 Numerical Results 

The whole set of algorithms presented in the previous sections has been translated into a 
FORM code and the output has been used to create a FORTRAN code. 

In this section we will present numerical results for two-loop two-point functions (self-energies), 
based on the results developed in the previous sections. In the following we present a detailed 
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comparison among different formulations for the same topology, different integration routines 
and, for few special cases, between our numerical results and the available analytical calculations. 
Typically we have used three different sub-routines for numerical integration: DOIEAF, a multi¬ 
dimensional adaptive quadrature over hyper-rectangle, multiple integrands; DOIGBF, a multi¬ 
dimensional quadrature over hyper-rectangle, Monte Carlo method; DOIGDF, a multi-dimensional 
quadrature, general product region, number-theoretic method. 


8.1 Numerical results for S 


121 


To check our results for this topology we consider two special cases, namely a) mi = m, ruj = 0 
for i >2 and b) m 2 = m^ = m,mi = m3 = 0 which correspond to Fq. (38) and to Fq.(45) of 
ref. [^. Results and comparisons for conhguration b) are shown in Tab. 2. The hrst entry is 
based on method 11, Fq. (0). and the NAG subroutine DOIFAF; the second entry refers to 


the same Eq.(TO) but uses the NAG subroutine DOIGBF. The latter usually returns a smaller, 
underestimated,integration error. The former has a larger error, typically 0.06% for the real part 
and 1.7%(0.9%, 0.5%) for the imaginary part at 81 GeV(83 GeV, 85 GeV). Relative deviations with 
the analytical results are < 0.006% for the real part and 0.09%(0.05%, 0.04%) for the imaginary 
part (0.07%, 0.02%, 0.01% for DOIGBF). As expected the largest errors occur around threshold. 

For conhguration a) we report results and comparisons in Tab. In Appendix E we give the 
analytical evaluation of the imaginary part of which receives a contribution from a two-particle 
cut and a three-particle cut. Gomparisons with the numerical result are shown in Tabs. The 
relative deviation with respect to the analytical result is < 0.04% over a wide range of energies, 
below and above the three-particle cut. In Tab. | we scan the threshold regions, y/s = m 3 + m^ 
and y/s = mi +m 2 + mi. Even at one MeV above threshold we are able to reproduce the imaginary 
part with a relative precision of less than 1% 


8.2 Numerical results for 


The interesting case for is represented by m3 = m^. No analytical results are usually 
shown in the literature since all authors write 


mi, m 2 , m 3 , m^, m 3 ) 


d 

dm3 


S'^^^(p^;mi,m2,m3,m4). 


(394) 


We have done the same and compared internally with the numerical derivative of 
However, numerical differentiation of functions that are already the result of numerical integration 
enjoy poor reputation, a fact that we conhrm. Unless a extremely accurate but time consuming 
evaluation of the multi-dimensional integral is performed we end up with large uncertainties in 
the derivative. Fortunately we can still check our results, at least for the imaginary part where 
an analytical approach is available. The results are shown in Tab. ^ and the larger error in the 
derivative, for larger energies, is due to the three-particle cut. The behavior of for m^ = m 3 
and around the normal threshold yfs = m 3 + m^ is shown in Fig. ^ and also reported in Tab. 
where we see that numerical precision is worst on the high energy side of the threshold, where 
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the imaginary part in non-zero. However, it is enough to go 100 MeV away from threshold to be 
better that 0.1% (0.3% at 50MeV). 


8.3 Numerical results for S 


221 


For we have compared our results with those of ref. il where the authors have introduced 
a one-dimensional integral representation of In particular we refer to Tab. 1 of ref. . 

Comparisons are shown in Tab. H where we have selected their numerical results with the smaller 
error. In this case we use the NAG subroutine DOIGDF which calculates an approximation to the 
integral using the Korobov-Gonroy number theoretic method with a Monte-Garlo error estimate; 
our two entries in Tab. ^ refer to the choice of 10193 ■ 101 and 22807 ■ 151 points respectively. 
Note that with their choice of masses, = i/nGeV, the normal thresholds are at 2.414 GeV, 
4.236 GeV (two-particle cuts) and 4.968 GeV, 5.146 GeV (three-particle cuts) respectively. Only 
the last energy is above all normal thresholds. There are two conhgurations where can be 
computed using the appropriate expansion: let variables x, y be defined as in Eq. (|421|) , we can give 
the \y\ -C 1 and the |x| -C 1 expansion of S'^^^(s; m, 0, 0, 0, 0), as discussed in Appendix F. Results 
are shown in Tab. By using Eq. ([424| ) we can compare the low-s behavior of m, m, m, 0, 0) 

which is shown in Tab. 0. 


9 Conclusions 


For one-loop diagrams the analytical way has been fully successful, resulting in compact ex¬ 
pressions containing at most di-logarithms For multi-loop diagrams, special cases have been equal 
successful, see for instance ref. and ref. [^, resulting in expressions containing generalized 
Nielsen poly-logarithms [^. However, this approach can hardly been extended beyond one-loop 
for general massive diagrams. For instance, the work of ref. and already shows that for 
two-loop integrals which contain a massive three-particle cut, the results in general cannot be 
expressed in terms of generalized poly-logarithms. Sometimes special cases of multi-loop diagrams 
can be expressed with the help of generalized hypergeometric functions or of other higher tran¬ 
scendental functions, but in the end this is equivalent to some series representation that always 
has a restricted region of convergence. In other words, from a practical point of view, they can 
hardly be extended above the closest normal or pseudo-threshold. 

This paper represents the second in a series devoted to numerical evaluation of multi-loop 
Feynman diagrams with the declared intent of covering all regions in the external parameters. In 
the first paper we have introduced the minimal Bernstein-Tkachov algorithm and discussed the 
simplest two-loop topology, or sunset diagram. 

In this paper we extended the previous work to cover all two-loop two-point functions, i.e. 
the self-energies. Each topology has been written with an integral representation in the Feynman 
parametric space which is particularly suited for numerical integration. For two external legs it 
turns out that distortion of the integration hyper-contour can be avoided. 
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Figure 21: The behavior of 5^^^, Fig. |^, around the normal threshold s = (m 3 + 
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Already dealing with the s un set topology we discovered this remarkable property that has 
been generalized now to all two-loop two-point functions but that remains true for all two-loop 
three-point planar topologies |^. For all these diagrams the matrix H of Eq.(|^) is singular. 
As a consequence, a change of variables is always possible in the quadratic form V so that the 
^-coefficient of Eq. (|I]) becomes independent on Feynman parameters, at least if we use only one 
iteration of the algorithm. In this way we are able to be closer to the original idea of the BT- 
method. However, as already stressed in the original paper ^j, this B will vanish at some threshold. 
We have analyzed all Landau singularities for the two-loop self-energies and have to show that this 
occurs at some non-leading Landau singularity of the diagram. Here some additional analytical 
work is needed before starting the numerical evaluation. 

A second motivation for this work was to start a comprehensive analysis of infrared divergent 
multi-loop diagrams, from the point of view of their numerical evaluation. Clearly the approach 
described in this paper is primarily intended for evaluation of massive multi-loop diagrams. How¬ 
ever, QED and QCD will be part of any realistic calculation and they usually lead to infrared 
singularities. Therefore, although the massless world is most efficiently treated in QED and QCD 
within the analytical approach, any multi-loop calculation in the standard model or beyond, being 
plagued by infrared divergences, requires to hnd a way of dealing with them in a purely numerical 
approach. 

For all two-loop two-point functions the relevant infrared objects are the on-shell derivative 
which are the building block of wave-function renormalization. For some of them the infrared 
divergent conhgurations are simple enough to allow for BT treatment with a consequent and 
straightforward analytical evaluation. The remaining ones, in particular when there is the pres¬ 
ence of overlapping infrared conhgurations, require a novel approach. Our solution is derived by 
adapting the general algorithm of ref. 0: the residues of the infrared poles and the hnite part 
of a multi-loop diagram can be cast into a form which allows for a reliable numerical integration. 

Work is in progress to extend the numerical evaluation to two-loop vertices, including non- 
infrared on-shell derivatives of self-energies, where the FORM codes for all topologies, including 
the non-planar one, have been developed . 
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10 Appendix A 


Consider the topology with mi = m 2 = 0 and m3 = m of Fig. |^. Let us consider the 
derivative 


TT 


4 ^ olll _ 
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(Bqi / crq2 


p-{q2+ p) 


ql (<?i - q 2 )‘^ (fe +pf + m 
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Figure 22: An S'^^^-like contribution to the W self-energy. 

We combine the first two propagators with a Feynman parameter x, perform the gi integration 
and combine the two resulting propagators with a new parameter y. After the g2-integration we 
put = —rri^ and obtain 

7“ l»= (^)' r (') B(1 - 1 - |) B(2 - 2 £, 1 + |), (396) 

where B is the Euler beta-function. The the on-shell derivative contains a single pole of ultraviolet 
origin and no infrared pole. 


11 Appendix B 

For completeness we report in this Appendix a list of master integrals diat are needed 
in the evaluation of see Eq.( p.21|) . Given X we introduce an auxiliary variable x defined as 
X = X — 1 and obtain 


dyy'^\ny = - 


dyy\Yi{l + Xy) = 
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dy{y + ln(l + Xy) = X t dy{l + Xy)-^ ln(l + Xy), 

Jo X Jo 

dy{l +Xy)-^ \n{l +Xy) = X, 

^ dy{y + ^)-Uny = U2 (^-^ dyy-Un{l + = - U2 
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12 Appendix C 

In this appendix we list all the integrals that are relevant to compute the on-shell behavior of 
sl^\ see Eq. (|227|) . We will introduce the following terminology: 

/*! r I —e/2 

X(/i, z/, a)= / dx x(l — x) J'{iJ,,p,a)= / dy y^ {1 + a y)'^. (398) 

Jo ^ ^ Jo 

We obtain the following results: 

77(--l,-l-e,X) = - 2X1(1-1-e, - ; I + 2’ 

= - x^’'"^ 2X1(1 -|-e, 1; 1+2’ 

The corresponding x integral gives 

J(--l,-l-e,X) = -^-^|i^ dxx'+^/2^1-x)"-^/2 
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2 r (1 + e/2) r (2 + e/2) “ T (n + 1 + e) T (n + 1 - e/2) 


r(l + e) 


E 

n=0 


(399) 
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A second class of integrals, 

[' dyy‘1^ (1 + A'»)-■-< = J(5, -1 - e, A), 
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are finite in the limit e —0. We obtain 
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A third class of integrals can be reduced to the previous ones: 


0 y + x 


= 1 . 


[ dx [ dy a;(l — x)] y"^"^ ^ {1 + Xy) = 

Jo Jo 

/ dx [ dy a;(l — |-(1 + X)“'^ + 2 lull + X) + e Li 2 (— X) — In^ 

JO JO L e L 

J dx — 2 lux + e Li 2 ^ + ln(l — x) lux | 


= -r 1 + - r 1 


where we used the relation 
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The constant is 


where we used 
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After integration we obtain 
where 
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(405) 

(406) 

(407) 

(408) 

(409) 
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13 Appendix D 


In this Appendix we report the explicit expression for the residne of the single infrared pole 
defined in Eq. (|305|) . 


221 


^ 221 =^ Jq dy{x2{l + xi)[{l + xif-{l-xiX2Yxly^]AQ^^ 


+ xly{l + xi) [(1 + xiY - (1 - xiX2y) y ] 

+ xi (1 + X 2 ) [(1 + X2f - (1 - xifxlxly"^] Aq2i 
+ xly{l + X2) [(1 + xa)^ - (1 - xiyfxjy'^] Aq22 
+ (1 + X2 - X 1 X 2 ) [(1 + X2 - xiXa)^ - xjxly"^] A^‘^ 

+ (1 + X1X2 - x^xa) [(1 + X1X2 - x^xa)^ - xjxjy"^] A^i 
+ xa (1 + xixay - XiXay) [(1 + xiXai/ - x^xai/)^ - xly'^] A^i + l} 


+ 


dxo 


dxi 

rl 


dy 


Xa rr(l-Xi)^ 5 i' 


1 + Xa 


{[ 


-2 


Xi 


+ ' 1 + Xa 


+ 


Bf-} 


poo pi 

- dx 2 dy —— 
2 Jo Jo 1 + Xa 


^(xi = 0) [xaLftdxi = 0) + 1] 


r*l 


POO p1 

+ / dxa / dxi I dy 

Ji Jo Jo 1 + Xa 


Xa cv{l — Xiy)^B. 


{[ 


-2 


+ 


(1 - xiyf 


1 + Xa 


bY 


1 

2 


dX2 


rl 


dxi 


n 


to 


1 + Xa 


B 2 {y = 0) [x2Lb2{y = 0) + 1 ] 


/ dx2 / dxi 
lo Jo Jo 


r 1 

r(l - XiX2yfB2il 

•-l + Xa 

y 


XiXa 


(1 + Xa)^ 

XiXa 

(1 + XiXa)^ 


(1 - xiX2y) B^i + 


Xa Al — xiy^B. 


22 ^- 


1 + XiXa L 


y 


2 Jo 


dX2 


{l-x,yfB^i} 

B2i{y = 0)' 


dxi {[■ 


(1 + X2)X2 -1 + 


1 + Xa 

3:2 


^21 ( 2 / = 0) [Lb2i{y = 0) + Inxa] + 


■ B22 iy — 0 ) 


L(1 + X1X2)X2 

■ 1 + X 1 X 2 ^ [Lh22{y = 0) + InXi]} 

1 A 

+ - dxi B^iiy = Xa = 0) [Lb2i{y = xa = 0) - 1] 

2 JO 


+ ^ / dx2 B^^{y = xi = 0) [Lb22{y = xi = 0) - 1] 

2 JO 


( 410 ) 


77 
































where the following functions have been introduced: 

^011 = (1 + Xif + (1 + xl){l + a;i)(l - xiX 2 )y + (1 - xiX 2 fxly‘^ 

^012 = (1 + xiYx2 + (1 + xly‘^){l + a;i)(l - xiX2y) + (1 - xiX2yfx2y^ 

Ao2i = (1 + X2Y + (1 + xlxl){l + X2){1 - xi)y + (1 - xiYxlxly‘^ 

A 022 = (1 + X2fxi + (1 + xlxly‘^){l + X2){1 - xiy) + (1 - xiyfxixly'^ 

Ai = {I + X2- xixlY + (1 + xl){l +X2- xix\)xiy + xjxjy'^ 

A 21 = (1 + X 1 X 2 - xlx2)^ + (1 + + X 1 X 2 - xlx2)y + xlxly^ 

A 22 = (1 + xiX2y - xlx2yfxl + (1 + xlxly‘^){l + xiX2y - xlx2y) + xlx2y‘^ 
Bi = 1 + X2 - X 1 X 2 + (1 - xifx2y 

82 = {I + X 2 - xiX 2 y)xi + (1 - xiyYx 2 

B 21 = {I + X 2 - xixly)xi + (1 - xiX2yf 

B 22 = (1 + X 1 X 2 - xlx2y) + (1 - xiyfx2 


and also 

1 13 

Laoii = In(Aoii) - - In y - ln(l + X2) + - ln(l - X1X2) - - ln(l + Xi) 

1 13 

La0i2 = ln(Aoi2) - -\nx2- ln(l + X2y) + - ln(l - XiX2y) - - ln(l + Xi) 

1 13 

La02i = ln(Ao2i) - -\ny- ln(l + X1X2) + - ln(l - Xi) - - ln(l + X2) 

1 13 

La 022 = ln(Ao 22 ) - - lu Xi - ln(l + XiX2y) + - ln(l - Xiy) - - ln(l + X2) 

11 3 

Lai = In(Ai) - - \ny + - Inxi + \nx2 - ln(l + X2) - - ln(l + X2 - Xixj) 

1 3 

La2i = ln(2l2i) - -\ny + Inxi - ln(l + X1X2) - - ln(l + X1X2 - x^X2) 

1 3 

J^a22 = ln(A22) + Inxi - - Inxs - ln(l + XiX2y) - - ln(l + XiX2y - x^X2y) 

11 3 

Lbi = ln( 5 i) + - Iny - - lna;2 + ln(l + X2) - - ln(l + X2 - X1X2) 

11 3 

Lb2 = ln(i?2) + - Inxi - - In0:2 + ln(l + X2) - - ln(l + X2 - XiX2y) 

1 3 

Lb2i = ln( 52 i) + - Inxi + ln(l + X2) - - ln(l + X2- Xixjy) 

1 3 

Lh22 ^ - - In X 2 + In(l + X1X2) - - ln(l + X1X2 - x{x2y). 


(411) 


(412) 
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14 Appendix E 


In this Appendix we consider all those results that are useful to derive analytical expressions 
for the imaginary parts of two-loop two-point functions. The two-particle cut contribution to the 
discontinuity of shown in Fig. 0, is given by 


Since the one-loop scalar function Bq is 


Bq = - 


where 


R = In 

pZ 


—m 3 ; mi, m 2 ) Imi?o (—s; m 3 , mi). 

(413) 

ln"^^7' R+^'^-^hn^l+2, 

jjL^ 2 p^ ml 

(414) 

— id + ml + ml — A 

i- - - -^-, 

2 mim 2 

(415) 


and where = A (— ml, ml), we must have the imaginary part of Bo (—s; m 3 , m 4 ) to O (e). It 
follows 


Bo (p 2 ; mi, m 2 ) = ^ - £ dx In ^ 

+ £{i(7 + lnx)Vic(2) + i/^‘ 


dx In 


A 

/i2 


7 -h In TT -h ^ In ^ 

2 /i^ 


}, (416) 


where, as usual, x = — ml)x + ml- Note that, for s > (m 3 -|- m 4 )^ and m^ > 

(mi - 1 - m 2 )^ the discontinuity develops a real part. However, in the same region the discontinuity 
coming from the three-particle cut of Fig. ^ also has a real part and the two cancel so that 
disc 2+3 5121 is purely imaginary. For completeness we report the result for disc 3 5^^^. Dehne hrst 

■s±± = (\/s ± "^4)^ - {m2 ± mi)‘^. 


2 _ 

r\j — - 


s++s__ 


and also 


Xs = —, Cl = ,/s++s__. 


— = 4 mim 2 s , — = 8 mim 2 (s — m, -f m? -|- m, -|- m\), 

Cl Cl 

C 5 A(s,mi,ml) 

— = 8 mim2 


(417) 


C4 (s - 

— = — 8 mim2 
Cl 


ms 


Cl 


ms 


It follows 


X4 (m2 — mi)^ X5 ml — (m2 — mi)^ 

X3 (m2 + mi)^’ X3 ml — (m2 + mi)^ 

2 in r 


(418) 


disc, 5^21 = 


Cl E(k) + C 2 K(k) + C3 n(a:3, k) 


+ C 4 n(a; 4 , k) -|- C 5 n(a ;5 — id, k) 9 (s — (mi -|- m 2 -l- m 4 )^), 


(419) 
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where E,K and 11 are elliptic integrals [^. The prescription is needed because in the last 
complete elliptic integral of third kind we have, in the region s > + and m\ > (mi+ 7712 )^, 

that xs > 1 and the integral must be understood as a Cauchy principal value. The corresponding 
imaginary part cancel the real part from disc 2 . Adding the two contributions we hnd 




2 L 


disc 2 + discs S 


121 


(420) 


For a numerical evaluation of elliptic integrals we use the NAG sub-routines S21BAF - S21BDF, 



Figure 23; The three-particle cut of diagram 
following the alternative dehnition by Carlson p5 |. 


15 Appendix F 

In this Appendix we report some special cases of for which the low-(high-)p^ expansion is 
known [^. Let us dehne the following variables 


m . ^ s . ^ 

x = -1-7 0 , y = - -—to. 

s 


We obtain the following two expansions: 

,2 


1 00 

S'2^^(s;m,0,0,0,0) = - [ 


n=l 


111 V + 2C(2) _ ^ ^ 


1 00 

.S2^^(s;m,0,0,0,0) = - y [ 


n=l 


where r 2 is dehned as 


2 n 

_ 4,,,+4+ 6 C( 3 )^ 

2n n'^ J s 


1 n-l 
= W 




1=1 


P 


Another known expansion is the following: 

OC 

m, m, m, 0, 0) = ^ 

2 m2 


n!F(3/2) / s W 

3 

In 7 / 

F(n-1-3/2) V4m2/ 

■{n + 1)2 

n + 1^ 


(421) 


|x| -C 1, 

(422) 

(423) 

(424) 
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yi [GeV] 

Re ^121 

ReS'^^^ (anal.) 

Im ^121 

(anal.) 

75 

54.27(3) 

54.2652(1) 

54.2624 

0.00 

0.00 

77 

53.11(3) 

53.1227(4) 

53.1076 

0.00 

0.00 

79 

51.63(3) 

51.6348(2) 

51.6333 

0.00 

0.00 

80 

50.59(3) 

50.6040(2) 

50.5868 

0.00 

0.00 

81 

49.35(3) 

49.3791(4) 

49.3415 

-0.58(1) 

-0.5850(5) 

-0.585457 

83 

47.68(3) 

47.6781(4) 

47.6763 

-2.18(2) 

-2.183(5) 

-2.18366 

85 

46.53(3) 

46.5088(4) 

46.5277 

-3.88(2) 

-3.880(5) 

-3.87974 


Table 1 : The topology for mi = m3 = 80 GeV and m2 = m^ = 0 as a function of ^/s. First entry is 
from Eg. ([7^ and NAG subroutine DOIEAF and the analytical result of ref. | 26 ]. Second entry is from 
Eq.(^) and NAG subroutine DOIBGF. The UV-pole is e = 1 and the unit of mass is /i = 1 GeV. 
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yi[GeV] 

Re ^121 

(anal.) 

Im ^121 

(anal.) 

75 

40.34(3) 

40.3501(4) 

40.3502 

-49.77(3) 

-49.7756(3) 

-49.7720 

77 

41.16(3) 

41.1952(5) 

41.1696 

-49.94(3) 

-49.9387(5) 

-49.9374 

79 

41.96(3) 

41.9512(4) 

41.9693 

-50.10(3) 

-50.1045(3) 

-50.0985 

80 

42.36(3) 

42.3447(4) 

42.3621 

-50.18(3) 

-50.1756(3) 

-50.1775 

81 

42.74(3) 

42.7302(4) 

42.7502 

-50.26(3) 

-50.2691(3) 

-50.2556 

83 

43.50(3) 

43.5399(4) 

43.5127 

-50.41(3) 

-50.4033(4) 

-50.4090 

85 

44.25(4) 

44.2784(4) 

44.2576 

-50.57(3) 

-50.5566(2) 

-50.5593 


Table 2: The topology for mi = 80 GeV and m 2 = m 3 = m^ = 0 as a function of ^/s. First entry is 
from Eg. ([ 7 ^ and NAG subroutine DOIEAF and the analytical result of ref. |26]. Second entry is from 
Eq.(^) and NAG subroutine DOIBGF. The UV-pole is e = 1 and the unit of mass is /i = 1 GeV. 
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[GeV] 



(anal.) 

above two-p cut 




180 

15.90(3) 

15.8909(5) 

15.8957 

190 

23.26(3) 

23.2562(5) 

23.2592 

200 

28.41(4) 

28.4148(4) 

28.4059 

210 

32.39(4) 

32.3947(5) 

32.3834 

220 

35.62(4) 

35.6124(4) 

35.6137 

230 

38.30(4) 

38.2994(4) 

38.2953 

240 

40.59(4) 

40.5769(5) 

40.5929 

250 

42.57(4) 

42.5749(5) 

42.5613 


Table 3: The imaginary part of for nii = m 3 = SOGeV and m 2 = m^ = 91 GeV as a function of ^/s. 
First entry is from Eq.(|^) and NAG subroutine DOIEAF. Second entry is from the same equation and 
NAG subroutine DOIBGF and third one is the analytical result of ref. |26| as reported in Appendix E. 
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[GeV] 


Im ^121 

(anal.) 

above three-p cut 




260 

44.30(5) 

44.3069(5) 

44.2915 

270 

45.84(4) 

45.8391(5) 

45.8339 

280 

47.24(4) 

47.2457(5) 

47.2316 


Table 4: The imaginary part of for mi = m 3 = 80 GeV and m 2 = m^i = 91 GeV as a function of ^/s. 
First entry is from Eq. dzi) and NAG subroutine DOIEAF. Second entry is from the same equation and 
NAG subroutine DOIBGF and third one is the analytical result of ref. |^| as reported in Appendix E. 
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a/s [GeV] 


Im 5121 (//) 

(anal.) 

I/ anal - 1 

Il/anal - 1 

around the two-p cut 






171.1 

1.26(1) 

1.25423(5) 

1.24533 

0.87% 

0.71 % 

171.2 

1.89(1) 

1.87835(5) 

1.87846 

0.13% 

0.01 % 

171.3 

2.38(1) 

2.37516(5) 

2.37086 

0.28 % 

0.18% 

171.4 

2.80(2) 

2.80222(5) 

2.79280 

0.43 % 

0.34 % 

around the three-p cut 






262.1 

44.63(4) 

44.6222(5) 

44.6337 

< 0.01 % 

0.03 % 

262.2 

44.65(4) 

44.6446(3) 

44.6392 

0.01 % 

0.03 % 

262.3 

44.67(4) 

44.6617(4) 

44.6589 

0.02 % 

0.01 % 

262.4 

44.68(4) 

44.6983(5) 

44.6718 

0.02 % 

0.06 % 


Table 5; A scan of the imaginary part of for mi = m 3 = 80 GeV and m 2 = m 4 = 91 GeV around the 
two- and three-particle cuts. First entry is from Eq.(79) and NAG subroutine DOIEAF. Second entry is 
from Eq.(|7^) and NAG subroutine DOIBGF and third one is the analytical result of ref. |^| as reported 
in Appendix E. 


85 
























[GeV] 

Im ^131 

0 

7 

-1.551(2) 

-1.551(3) 

12 

-0.3047(3) 

-0.3047(4) 

17 

-0.1494(1) 

-0.1494(11) 

22 

-0.09166(9) 

-0.09183(282) 


Table 6 ; A comparison between 
mi = 1712 = ^4 = 5 GeV as a function 


and the numerical derivative of for m 3 = ms = 
of -v/i. The UV-pole is e = 1 and the unit of mass is ^ 


1 GeV and 
= 1 GeV. 




Re ^131 im ^131 


d Tv^ Cl21 


2.46581(1) 

0.00 

0.00 

2.76299(1) 

0.00 

0.00 

3.20714(1) 

0.00 

0.00 

3.99512(2) 

0.00 

0.00 

3.21(1) 

-5.28(1) 

-5.2794(5) 

1.820(2) 

-4.381(4) 

-4.3735(8) 

1 . 202 ( 2 ) 

-3.834(1) 

-3.8354(3) 

0.8387(5) 

-3.458(1) 

-3.4588(5) 

0.5954(5) 

-3.1725(7) 

-3.1719(4) 


Table 7; and the numerical derivative of for m 3 = ms = 1 GeV and mi = m 2 = 
around the normal threshold ^/s = 6 GeV. The UV-pole is e = 1 and the unit of mass is fj, 







s [GeV^] 

Re ^221 

ReS^^hel 27 

Im ^221 

Im^221^gf_ 27 

0.1 

-0.287205±0.20 x lO''^ 

-0.287237±0.35 x 10“® 

-0.287238±0.25 x 10“® 

0.00 

0.00 

0.5 

-0.294578±0.15 x lO""* 

-0.294590±0.23 x 10“® 

-0.294952±0.25 x 10“® 

0.00 

0.00 

1.0 

-0.304543±0.21 x lO''^ 

-0.304522±0.22 X 10“® 

-0.304521±0.25 x 10“® 

0.00 

0.00 

5.0 

-0.452499±0.17 x lO''^ 

-0.452518±0.31 x lO”"* 

-0.452520±0.28 x 10“® 

0.00 

0.00 

10.0 

-0.488036±0.59 x lO"® 

-0.488224±0.18 x 10“® 

-0.488153±0.23 x 10"® 

-0.353197±0.69 x 10“® 

-0.353082±0.18 x 10“® 

-0.353217±0.23 x 10-® 

50.0 

0.173897±0.49 x IQ-^ 

0.173900±0.18 X 10-4 

0.173901±0.22 X 10“® 

-0.118089±0.40 X 10-4 

-0.118082±0.20 X 10-4 

-0.118080±0.22 X 10-® 


Table 8: compared with the numerical result of Tab. 1 in ref. [27|. Here = nGeV^, for 

n = 1,...,5. The second entry refers to a run with a much larger number of points in sub-routine 
DOIGDF. 




























s [GeV^] 

Re 5221 

Re S'221 expansion 

Ini^22i 

Im 5*22^ expansion 

0.1 

0.3740(6) 

0.373988 

0.27957(9) 

0.279970 

1.0 

0.19599(4) 

0.195930 

0.20847(6) 

0.208394 

10.0 

0.07144(3) 

0.071467 

0.14067(3) 

0.140692 

5.0 

-1.61744(9) 

-1.61746 

0.4909(2) 

0.490807 

10.0 

-0.79270(2) 

-0.79270 

0.14061(5) 

0.140692 

50.0 

-0.150366(2) 

-0.150362 

0.007492(3) 

0.007492 


Table 9: compared with the expansions described in Appendix F for mi = 100GeV(mi = 1 GeV) 

upper (lower) part of the table. The sum over n in Eq. (|422D is restricted to n < 20. 

















s [GeV^] 

Re 5221 

Re S'221 expansion 

Ini522i 

Im 5*22^ expansion 

0.1 

0.0496(2) 

0.049542 

0.01566(3) 

0.015709 

0.5 

0.04154(5) 

0.041506 

0.015722(8) 

0.0157145 

1.0 

0.03805(3) 

0.038051 

0.015717(9) 

0.0157211 

5.0 

0.030085(9) 

0.0300729 

0.015776(4) 

0.0157739 

10.0 

0.026682(9) 

0.0266732 

0.015836(6) 

0.0158406 


Table 10: compared with the expansions described in Appendix F for mi = m 2 = m 3 = 100 GeV 

and m 4 = ms = 0. The sum over n in Eq.(424) is restricted to n < 20. 
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